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ABSTRACT 


There has been a reneuied interest in the time-dependant 
quantum mechanical (TDQM) approach to reactive scattering in 
the last feui years. The different computational methods and 
applications that have created such an interest have been 
reviewed in Chapter 1 of the thesis. A comparative study of 
the efficacy of 3-point, 5-point and 7-point finite difference 
formulae and f ast- fourier-transfo rm algorithm to compute the 
Laplacian operator in the Hamiltonian of the Schrodinger equa- 
tion has been made and the results are reported. 

The absorption spectrum of the transition state (TS) con- 
figurations for the collinear reaction 

H + ^ + H (R1) 

has been predicted by solving the time-dependent Schrodinger 
equation for the nuclear motion on the chemically accurate 
Siegbahn-Liu-Truhlar-Horou/itz (SLTH) potential-energy surface 
(PES) in Chapter 2. The u/ings to the Lyman-d line have been 
interpreted in terms of the build up of the probability density 
function in the inner repulsive u/all of the PES, Effects of 
varying relative translational energy of the reactants and the 
vibrational state of the molecule under attack have been invest 
igated. Transition state spectra have also been computed for 
thermal distributions of H + at 300 and 1000 K, Our result 

are in agreement with those obtained by using the tirae-independ 
ent quantum mechanical approach as u/ell as the quasiclassic al 
trajectory||(QCT) method. 



Dynamics on the PES that correlates with H (2p) + H 2 
has also been carried out using the TDQM approach and the results 
are reported in Chapter 3. The emission spectrum of the TS con- 
figurations for the collinear reaction (Rl) has also been comput- 
ed and it is shown to be complementary to the absorption spectrum 
results reported in Chapter 2. 

Results of a detailed inv/estigation of the collision- 
induced dissociation process in the collinear H + collisions 
on the SLTH surface, using the TDQM approach are given in Chapter 4 
The exchange and the dissociation probabilities have been computed 
over a wide range of relative translational energies (0.5 -10,0 oU ) 
The dynamical threshold for dissociation ('^8-10 eV) has been found 
to be nearly two times the energetic threshold 4,5 eU), This 
is in contrast to the results of three dimensional QCT calcula- 
tions carried out elsewhere (3. Dove et al.) for the same system 
wherein the dynamical threshold was comparable to the energetic 
threshold. But, our results compare well with the earlier reported 
TDQM results for collinear collisions on a model PES, Therefore, 
we have concluded that the difference in the dynamical and energe- 
tic threshold arises mainly from the constraint that the three 
nuclei lie along a straight line. For v * 1 of H 2 molecule, the 
nodal structure in the initial wave function has bsen found to bo 
retained, even after the collision was over. In addition, more 
nodal structures in the wavepacket densities became apparent 
suggesting that this collisional process is governed by a propen- 
sity rule of Av = 2 for vibrational excitation, 

A summary of our findings is presented in Chapter 5* 
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CHAPTER 1 


INTRODUCTION 

1 . 1 I ntroduction 

Although it u;as realised more than fifty years ago that, 

in principle, any chemical problem could be solved by an appro- 

-] 

priate formulation of it in quantum mechanical terms, in 

practice, developments in that direction had to au/ait the 

arrival of modern electronic computers. In the early seventies 

R 

one could calculate reaction probabilities (P ) readily for 
exchange reactions of the type 

A + BC > AB + C (R.1) 

in collinear geometries by solving the time-independent 
Schrddinger equation for the nuclear motion, given the elec- 
tronic energy (for example, see ref, 2), That exact (converged) 



2 


solutions to the problem in three dimensions (3D) also could be 

3a b 

obtained tuas demonstrated in 1975 mhen Kuppermann and Schatz ' 

3c 

and Elkoiuitz and Wyatt published their results for the 
reaction 

H + H 2 ^ H 2 + H ■ (R2) 

at louj energies, making full use of the symmetry of the system. 

There has not been any report of results of similar accuracy for 

any other reactive system in the past ten years. Even for the 

4 

reaction (R2) , some converged results have become available 

only recently on the chemically accurate Siegbahn-Liu-Truhlar- 
5 

Horouiitz (SLTH) potential-energy surface (PES) for the ground 
and the first excited vibrational states (v = 0, 1) of the 
reactant. While some progress can be expected in the next feui 
years, thanks to the supercomputers/parallel processors becom- 
ing available, the enormity of the problem remains, 

A viable alternative is a semiclassical (SC) approach. 

"For any system u/ith a classical analogue, a state for which 
the classical description is valid as an approximation, is 
represented in quantum mechanics by a wave packet. The wave 
packet should move according to the laws of classical dynamics'* 
(P.A.M. Dirac). ^ Equations of motion were derived by Dirac, 
assuming "h to be very small as in the limit "fi 0 , the classi- 
cal dynamics corresponds to the quantum dynamics, Heller has 
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taken a slightly different v/ieuipoint in the series of papers 

he has published. ” He has exploited the correspondence 

betiyeen the classical dynamics and the quantum dynamics of 

localised ujave packets (WP) on smooth PESs, ' He assumed that 

if the luavef unction is represented by a sum of Gaussians, u/e 

can propagate each Gaussian independently and that in the process. 

the width of each Gaussian is smaller than the length over which 

the potential changes. He has applied successfully his semi- 

0-11 


classical method in the study of photodissociation. 


photo- 


12 13 1 

absorption, molecular spectroscopy, stochastic processes, 

15 1 6 

investigation of local modes and dynamical tunneling * and 

17 1 8 

atom diffraction by surfaces, * His studies have also sti- 

19 20 

mulated elegant experiments, * 

21 

Coalson and Karplus have extended the Heller’s method 

(HM) to provide a theory of quantum mechanical (QM) motion for 

multidimensional anharmonic systems. It has also been extend- 
22 

ed for systems for which the generalised Liouville equation 

is applicable. Its description of the atom diffraction from 

1 7 1 8 2 2 b 

surfaces ’ has recently been improved by Oackson and Matiu 

in a way that would enhance the prospects of its contribution 

to surface science. 


23 24 

It has been shown recently * that the HM is valid 

only for systems which deal with short-time dynamics of quantal 

degrees of freedom and for high-energy processes. Recent cal- 

25 26 

eolations by Skodje and Truhlar^ and Heather et al, indicate 
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that the HM fails to giue correct values for the tlnie evolu- 
tion of the states of the Morse oscillator. They have shoujn 
that because of the anharmonicity of the oscillator potential 
in the asymptotic region, error in the time evolution accumu- 
lates and this causes large inaccuracies in the transition 
probabilities especially for lou/-energy collisions, ufhich 
require a long time to establish the asymptotio conditions. 

For many elementary systems, it has been possible for 

sometime to solve the time-dependent SchrSdinger equation 

(TDSE) directly (numerically) by time-evolving a WP and 

thereby obtain a QM description analogous to the classical 

27 

picture. In their pioneering work, Mazur and Rubin solved 

the TDSE for a reaction of the type (R1) using a model poten- 

2 8 

tial based on Harmuth’s integration scheme, McCullough and 
29 

Wyatt (MW) did an elaborate study in luhich they solved the 

TDSE for the reaction (R2) in collinear geometries using a 

realistic PES,^*^ Unfortunately their approach could not be 

used readily for many systems as (i) the computer memory and 

the time requirement were large, (ii) the resolution of the 

stat e-to-state P information ujas poor and (iii) there mas 

very little hope of the method being extended to 3D, In the 

last feuj years, houiever, the situation has changed dramatic- 

31 32 

ally, Neujer methods * of solving the TDSE have become 

available and they are an order of magnitude faster than that 
29 3 3 

used by MW. Felt ot al, have proposed another method of 
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solving the TDSE, based on the spectral properties of it, for 

cases uihich deal with bound state problems. Thus, there is a 

renewed hope of studying 3D atom-diatom exchange reactions 

34 35 

quantum mechanically. Also, Kulander and Leforestier have 
shown that it is possible to deconvolute the average reaction 
probabilities ( <P > ) that result from a time-dependent quantum 
mechanical (TDQM) study. They have also shown that the TDQM 
approach can be used successfully in solving collision-induced 
dissociation (CID) processes - an area where the time-indepen- 
dent approach has been carried out successfully only for a few 
model systems. There are also other advantages offered by 
the time-dependent approach. It is analogous to classical 
trajectories in that a pictorical representation of the chemical 

reaction is possible-giving us a "feel" for the molecular colli- 

29 

sion. It has been shown that flux patterns could be obtained 
for a chemical reaction on a PES and the whirlpools in the 
resulting picture can provide insight into reactive scattering 
resonances. Very recently, we have also shown that by solv- 
ing the TDSE, we can construct the probability density function 
and predict the absorption spectrum for the transition state 
(TS) of the collinear reaction (R2) , We have also predicted 
the emission spectrum for the TS of the same reaction on an 

expited state PES, Gas-surface scattering problems have also 
38 

been shown to be amenable to, solution via the TDSE. 
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The TDQm methods haue been applied to different oollinear 
triatomic systems. The only exception is the study of Dackson 
and Wyatt who studied the time ev/olution of a planar 
system. The speed of the explicit method gives us the hope of 
extending it to atom-diatom exchange reactions in three- dimen - 
sions in the very near future. The existing methodology and 
applications of the TDQM approach to a study of reactive 
scattering processes are reviewed below. 


1.2 Methodology 

97 

Mazur and Rubin solved the TDSE 


i'Ka^ ( r jt)/ a t = Hf (^,t) .. ('I) 

2 

where H = - ^ + U ( r) (2) 

2m ' 

with ^ the Laplacian operator and \l(r) the potential energy, 
by separating the real and imaginary parts in 'll . Letting 
^ = u + iv, Eq. (1) becomes 

Hu = --h ^ 

nu -n 3 ^ 

Hv = 'K .. (3) 

The spatial and the temporal derivatives were approximat- 
ed by second order and first order finite differences (FD) , 
respectively. 
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With the advent of high speed computers, later workers 
could deal with the complex wavefunetion directly. The time 
evolution of the wavefunetion can be represented as 

^(t) .. (4) 

with U = exp [ -iAtH/ti] , At = t - t^ .. (5) 

li(t) is the wavefunetion at any desired time t and^(t^) is 

the initial wavefunetion at time t^. For computational reasons, 
the (exponential) evolution operator in Eq, (5) has to be trun- 
cated, but at the same time the unitarity of the operator has 

29 

to be preserved. Therefore, MW represented it as 

U = [1 -iAtH/ii] /[1 ■••iAtH/ln] .. (6) 

When time evolution is carried out in steps, ff at the n and 
n + 1 steps are related as 

= U'i'"' •• (7) 

and .. (8) 

Eliminating between Eq, (7) and (8) gives 




.. ( 9 ) 
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Thus, 

[l-iAtH/'h] [ 1+iAtH/in] 

Employing a 3-point ED scheme for the second order spatial 
derivative of ^ in one dimension gives: 




n+1 4-2 n+1 n^1 n+1 - n+l 


UihAre j refers to the value of ’J? at the location x^, that is, 


ifj'’ =S(Xj, t,) 

T 

Substituting for in Eg, (10), uje get 


.. ( 12 ) 


n^1 


n+1 n+1 


[1 + i(2a + \/jAt)] ^j-1 


) = 


[l-i(2a + UjAt)] 


where a = 'KAt/2m(Ax) 


.. (13) 
.. (14) 


Defining “ as a vector with components 5?. , Eq, (13) 

0 


becomes 


n+1 n«1 

kll = 8W 


.. (15) 


where = A and the elements of ^ and _B arp easily deduced 
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from Eq. ('13), The above method is referred to as an implicit 
method and is computationally expensive. It requires a large 

computer memory as ivell. As a result, it has been applied to 

r. 4. n 40_42 

a feuj systems only. 

31 

Askar and Cakmak have proposed an integration scheme 

for solving the TD3E in ujhich U is truncated after the linear 

term. By combining the forward and backward evolution opera- 

tors, they could devise an algorithm which is an order of 

2 9 

magnitude faster than the implicit method of MW, 

Subtracting (8) from (7), we get 




[ exp{ -iAtn/h } - exp {iAtH/h}]^*^ 


.. ( 16 ) 


Expanding exp[+iAtH/li3 into a Taylor series and keeping only 
the linear terms, we obtain 




n+1 


= - 




n-1 


.. ( 17 ) 


The second order spatial derivative of W for an equally 
spaced mesh in x, using the 3-point FD formula, is 


dx2 


n n 

-J- :!. 


2 * 


Ax' 


. . ( 18 ) 
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Substituting Eq. (18) into Eq, (17), uie get 


n+ 1 n n n 1 

= - 2 i[( 2 a + \/.AtA)fj^-a(’ 1 '^^^ .. (I 9 ) 

This can be ujritten as 


n+1 




( 20 ) 


This method is referred to as an explicit method, as the state 
of the system at any time step (n+1) is directly calculated in 
terms of the states at past times without the need to invert 
a matrix. 

In two dimensions, Eq. (I 9 ) becomes 


n+1 n- 1 


^j,k ' . " 2a.[{ 2(a-, +a2) +Uj^^AtA} 

J f 




)] 


.. ( 21 ) 


It turns out that the above first order explicit method 
is identical to a second order explicit method* Therefore, the 
errors are of third order in At, As a matter of fact it can 
be shown readily that the turncation upto any odd order in 
time is equivalent to the truncation upto the next higher 
order. 
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Truncating the foriuard evolution operator upto the 
second order term gives 

U(At) = [l -i&tHAi + (iAtn/H )^/2!] (22) 

Truncating the backiuard evolution operator upto the second 
order term gives 

U(At) = [l +iAtH/'h + (iAtH/li)V2! ] .. (23) 

Substituting for the forward and backward evolution operators 
from Eqs. (22) and (23) in Eq. ("IS), we obtain 

^n+1 ^n-1 ^ [.2iAtH/li] (24) 

u/hich is identical to (17), 

34 35 

TDQM studies ’ have also been carried out using 
either a 5-point FD or a 7-point FD scheme instead of the 
3-point FD scheme [Eq,(l8)] to compute the Laplacian operator 
in the Hamiltonian. We have carried out a comparative study 
of the accuracy of the 3, 5 and 7-point FD schemes for evalu- 
ating the second derivative of a function using a Gaussian, 
expf-(x-x^) "] as a tost case. The rms deviation is plotted 
as a function of the mesh size Ax, in Fig, 1, It is clear 
that the 3-point FD can yield reliable results if Ax is 


RMS DEVIATION (xl6^) 


V - 3 POINT FD 
X - 5 POINT FD 
• - 7 POINT FD 


0.10 


0.20 


Pig. 1. RMS deviation of the second derivative of the function 

2 

exp[-(x-x^) j plotted as a function of the mesh size 
for different FD schemes. In all the cases, 0^x<_6.3 
with X = 3 . 15 
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sufficiently small. On the other hand, the higher order FD 
methods allouj the use of successively largerAx for obtaining 
results of the same accuracy. Results of our comparative 
study of these methods for the reaction (R2) are presented in 
Chapter 2, 

32 

Recently, Kosloff and Kosloff have proposed a Fourier 

transform (FT) method to obtain the kinetic energy term. They 

make use of the property of the Fourier transform that a 

derivative in the spatial domain becomes a multiplication by 

i[<ij^ in the frequency domain, inhere is the wave number, 

corresponding to the 1-spatial coordinate. The Laplacian in 

Eq, (2) is thus obtained by performing a 2-dimBnsional Fourier 

transform on 5? , multiplying the result by -(K. + K. ) and 

^1 ^2 

performing an inverse Fourier transform back to the spatial 
domain. In practice, the readily available, fast Fourier 
transform (FFT) algorithms are made use of. 

32 

Kosloff and Kosloff approximated the temporal deriva- 
tive of by the second order FD scheme 

3^*^ ^ n+ 1 ^n— 1 

_ = 2 At 

3 1 

which is equivalent to the first order explicit method for 
the time evolution. 
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Substituting Eq. (25) in Eq, (1), uie obtain 


sn-^1 _ 5,n-1 
2it 


.. ( 26 ) 


Rearranging the above equation, uie get 
^n+1^ . 2iAtHlg^ ^ ^n-1 

which is the same as Eq. (17). 

32 31 

Therefore, the FT method and the explicit method 

differ from each other only in the way in which the kinetic 

energy term is computed. The disadvantage of using the FD 

method is its low accuracy and slow rate of convergence, due 

to the semilocal approximation to the kinetic energy. Because 

the FT method uses a nonlocal function expansion, the number 

of operations per grid point is higher than for the FD method, 

43 

But, by using the FFT algorithm, the numerical efficiency 
can be considerably improved. Since the FT method is completely 
vectorisable, its efficiency can be improved using array 
processors. 

We have also made use of the FT method to solve the TDSE 
for the reaction (^2) , the results of which are reported in 
Chapter 2. We have checked our computer code 

by computing second derivatives of the function Gxp[-(x2+y^)] 
where -5_<x<5 and -5< y<5, with a mesh size of Ax = Ay = 0 . 125 . 


.. (27) 
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-4 

Ws obtained an rms deviation of 3.8 xlO in contrast to the 
value of 2,0 X 10”^ obtained by the 7-point FO scheme. There- 
fore, uie see that the FT method is an order of magnitude more 
accurate than the FD scheme for the same mesh size. The 
function and its Laplacian are shown in Figs. 2 and 3, respect- 
ively. Because of the greater accuracy of the FT method, we 
can afford to use a larger mesh size, when compared to the FD 
scheme. This would result in considerable reduction of the 
number of grid points to be used for dach spatial dimension. 

The above fact would play an important role in carrying out 
3D quantal calculations. 

The FT method has generated a lot of interest in the 

44 -. 49 

time-dependent approach to reactive scattering. Tal-Fzer 

45 

and Kosloff have proposed a new propagator scheme for solv- 
ing the TOSE in which the evolution operator has been represent 
ed in terms of a Chebyschov series: 

N 

BXp[-iAtH/hj = S a^$^ (-iAtH/+i) .. (28) 

n=0 

where a„ are the expansion coefficients and § are the 

complex Chebyschev polynomials. This new method for the time 

evolution combined with the FT method for the evaluation of 

the Laplacian seems to have an error determined only by hard- 

45 

ware considerations. 
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( 5 , 5 ) 



(- 5 ,- 5 ) 

Fig. 2. Plot of the test function exp[- (x^+y^) ] with 

”5^y^5 
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( 5 , 5 ) 



Pig. 3. Laplacian of the test function 
plotted in Fig. 2 ■ 


Recently, Bisseling and Kosloff"^^ have proposed yet 

another method for the solution of the TDSE expressed in 

spherical or polar coordinates. The radial part of the 

Laplacian has been computed using a Fast Hankel Transform 

43 

method , based on the FFT algorithm. It has been success- 
fully applied^^ for the collinear reactive system in Delves 
hyper spherical coordinates. 

In order to compute the eigenvalues and eigenfunctions 

33 

of the Schrodinger equation, Feit et al. have proposed a 
method, based on the spectral properties of the TDSE. The 
method requires the computation of a correlation function 
(r,t^)llif (r,t)> from a numerical solution off(r,t). They 
have adopted a different strategy to compute SC'fjt), given 
the initial luavef unction r,t^) . 

Since 

f (r,t^ + At) u(At)5'rr,t^) .. (29) 

and u(At) = exp -i^] ..(3D) 


^(r,t + 

' » o 




.. (31) 


Similarity transformation of U leads to. 

At) = S-'' exp[i^i||^].axp[-i^] S*(*r,t^) .. (32) 



Substitution of 


s = axp[.iA|E£], S--' = ..(33) 


in Eq, (32) results in 


!E(?,t„+At) = pxp[i^^] sxp[-i^] 


IfJi] .xp[i^||l^].5^(?,t„) 


0-^ 


.. (34) 

33 

Equation (34) forms the split operator algorithm for advanc- 
ing the ujauef unction by an incremental time At» 

The operator, exp [ (■^‘^^) j applied to^Crjt^) is 
equivalent to solving the free particle luave equation 


at 




.. (35) 


over a time At/2 uiith ( r ,t„) as the initial uiavef unction 
at t = t^. The solution to Eq. (34) is obtained ujith the 
help of the band-limited Fourier series representation 


^(r,t) = 


n/2 

E 

-N 


N/2 

S 

-N 


n/2 

^ ^jkl exp [t^( j" +l<y ♦ Izl] 


j =-^+1 |<=-^1 



20 


uihere is the length of a side of the computational grid. 

The RHS of Eq, (34) is thus equivalent to free particle propa- 
gation over a half-time increment, a phase change from the 
action of the potential applied over the u/hole time increment, 
and an additional free particle propagation over a half-time 
increment . 

This procedure has been shown to be accurate and effi- 

43 

cient when implemented with the help of an FFT algorithm. The 

sampling interval At must be small enough to accommodate 

atleast the entire spectrum of bound state energy levels for 

the system. The effectiveness of this spectral method has 

33 

been demonstrated in the case of a one-dimensional double 

well potential and for the two dimensional Henon-Heiles poten- 

„ 50 

tial* Subsequently its utility has also been tested for 

practical three-dimensional potentials applicable to 502, 0^ 

and H 2 O, It has also been used to time evolve a Wp in the 

t 

Henon-Heiles potential in order to identify chaotic behaviour 
51 

of the system, 

1,2,1 Stability criteria 

A good integrating scheme to solve the TDSE has to be 
stable to round-off errors in the computation. Yet it has to 
be general enough to deal with most of the physical potentials 
and initial wavepackets. Also, the numerical algorithm has 
to be reasonably fast. The stability criteria to be met by 
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the choice of steps in time and space for the different methods 

2 8 31 52 53 

discussed above have been worked out elsewhere^ ’ » ♦ Only 


the essential features are outlined below. 

If the error eP = exp(iqx-) in is magnified by 
times in the next step, that is, e. = 9 > then for the 

method to be stable, Igl should not exceed unity. The ratio 
g = e / is called the growth factor. 


n+1 


Substituting 'P. + and W + (g)— eP for 

J J J J 


and 


in Eq. (19), we obtain 


g + 2i[2aCl-cos( qAx)}+ V -At/iS] g - 1 =0 

3 


.. (37) 


where a is defined as in Eq. (1A), 

The magnitude of each of the two roots of the Eq. (37) 
would be equal to 1 if the following condition is satisfied; 


1 > {2a[ 1-cos( qAx)] + V .At/in J 

vJ 


.. (38) 


As q may have any value, Eq» (38) leads to the following condi- 
tions for the stability of the integration; 


1 i 4a ? At< (1-4a)+i/|V . I 

“ ' J max 


.. ( 39 ) 
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In tuio dimensions, 

1 > 4 (a^ +a2) j At< [1 -4(a^ .. (40) 

The TDSE resembles the heat diffusion equation in that 

it involves a first order temporal derLvative and a second 

order spatial derivative. Yet the former seems more amenable 

to numerical solution than the latter. The imaginary number 

i in the TDSE changes the character of the differential equation 

3 1 

and results in the corresponding numerical scheme being stable. 

The stability conditions of the FF method have been uuorked 

32 

out for the case of a free particle system. Houjever, with a 
knowledge of the physics of the problem, one can come up with 
optimum values of step size in time and space for cases which 
deal with either a potential barrier or a potential well, 

1,2,2 Choice of initial conditions 

To carry out the dynamics of an atom-diatom exchange 

reaction of the type (R1) , if wie restrict outselves to collinear 

collisions, only two coordinates r-^ (the A-B distance) and r 2 

(B-C distance) are needed to specify the geometry. It is 

54 , 

possible to choose a set of coordinates (q-j, q 2 ) vi^hich reduces 
the motion of three particles in ( , r^) coordinates to the 
motion of a single particle, nost of th»e studies so far have 
omoloyed the scaled and skewed coordinates (q>|, ^ 2 ) which 
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diagonalize the Hamiltonian 


H = -(4iV2/J^^g2)C3VaqT^ + d^/dq^^) +\y(q^,q2) .. (41) 

A linear transformation that effects the required diagnoliza- 
tion is 


q^ = r^ 0 T 2 sin© 
q 2 = /3r2 cos© 


u;ith 

and 


0 = [i mgCm^ + mg)} /{m^(mg + m^) }] 


1/2 


0= sin” [m^mg/l(m^ + mg) (mg + mg)}] 


1/2 


.. (42) 
.. (43) 

.. (44) 

.. (45) 


where m^ , mg, mg denote the masses of the atoms A, B and C, 
respectively. 


The coordinate system is shown schematically in Fig, 4. 
In such a coordinate system, ^(t = 0) is chosen as 


*(t = 0) =S(q-,,q2) = *(qi) .Xn('l 2 ) 


.. (46) 


Where is the n-th state vibrational eigenfunction of the 
diatom BC; 

[-('1^V2AJgg) 3V3q2^ + V(q^ = qd , q2)] X^(q2) = E^X^(q2) 


. (^ 7 ) 
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For a Morse oscillator, the ground state eigenfunction is 
giv/en by 

.. (48) 

u/here is the normalization constant for the ground vibra- 
tional state given below, b = [( 8i;D)^af,].1 , u. is the reduced 
mass of 8C, D is the dissociation energy of 0C, a is the Morse 
parameter and is the q 2 -space value for the equilibrium 

internuclear distance of BC. The normalization constant for the 
n-th vibrational state Morse wavef unction^^ is given by 

= ..(49) 

2 

luhere M Q ^ S ■^(A-2n‘*-S-1) C 9 n\ 

n „«(A-2n-1) SI -• (50) 

where A = b^l and r represents the gamma function. 

The first vibrational state wavefunction for a Morse 
oscillator is written as 

X, (qj) = N^F^FjF, ..( 51 , 

inhere Fi = exp(-DF)i Fj = f((*“^'''^5'^2) ^ 

Fj = Z".(A-n-1)nz(''-''); Z = AF , D = a/2 , 

F = exp [ -a(qj . q 2 °)J; 
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and is the normalisation constant for the first vibrational 
state. The first vibrational state uuavefunction has a node 
located at the point ^2 ~ '^ 2 °* 

The question of the validity of representing the UJave- 
function uiithout any concern for the relative phase of the 

vibrational eigenfunction for the time-dependent Wp calcula- 

27 29 53 56 

tions * * has been recently addressed by Agraitial et al. 

They computed the rate coefficient for the collinear reaction 


T + HD ^ TH + D 


(R3) 


53 

using the explicit method, ujith and luithout the phase factors 
of the vibrational eigenfunction being incorporated. The results 
mere identical in both cases, thus providing justification of 
the earlier approach. 

$(q-j) in Eq, (41) is taken to be a WP constructed by 
superimposing many planeivaves to represent the translational 
motion of the atom A with respect to the center of mass of BC, 
Typically, it is taken to be a minimum uncertainty Gaussian Wp 
in the momentum (k) space, centered at k^^ corresponding to the 
classical relative translational energy C^tr> of the reactants 
such that 


^tr ='h^C''2‘'A,BC 


*( 111 ) .. (52) 
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uihere is the center of the Wp in the coordinate space. An 

optimal width (6) of the Wp has to be employed in the computa- 
tions, In principle, if one chooses a Wp whose k space distri- 
bution approaches a delta function, then the results can be 
directly compared to the fixed collision energy results obtained 
either via QCT or time-independent methods. In practice, how- 
ever, Such a procedure is not feasible since the resulting Wp 
will be completely diffused in the coordinate space. This would 
result in a non-physical situation, in which even at t =0, the 
reactants would not be well defined. Hence, the integration of 
the TDSE has to be carried out over an extensive configuration 
space grid such that the atom A is separated well from the mole- 
cule BC at t=0 requiring a substantial amount of computer time. 
On the other hand, if the Wp has a narrow distribution in the 

configuration space, it becomes diffused in the k space, thus 

R 

rendering any determination of the energy dependence of P 

difficult. In most of the studies so far the value of 6 has 

-1 

been chosen to be 0,25 a,u. 

If we are interested in a state (v)-selBctBd rate coeffi- 
cient at a translational temperature T, the following Wp has to 
be chosen; 

$(qi) * [t(3A) iF^£ 3/4, 1/2, -0^(q^-q° >2} 

.. ( 53 ) 

--2ir(5/4)(q^-q^°)0 /i^5/4, 3 / 2 , -3^(q^ 
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Uihere represents the confluent hypergeometric function and 

/ 2 1/2 

3 = (tfkT/1i ) . The of the above Wp corresponds to 

i' 2 2 

211 0 /li. This Wp had been previously utilized by Mazur and 
27 29 

Rubin and MW, If an overall rate constant is desired, an 
appropriate summation over the v states can also be incorporated 
in the initial ujavef unction, in the Eq. (46). 

For computational purposes Eq, (46) has to be discretised 
and the values of the initial mavef unction have to b e generated 
over a set of points in the - plane given by q>, = (j-'i)Aq^ 

+ t ^2 ~ (i-'I)Aq 2 ^2° ^here j = 1, 2, ..,,3, and i = 

1, 2,, ..I, called a grid and q^° and q^^ are convenient starting 
values for q/| and q 2 , respectively; Aq-| and Aq 2 are small 
increments and determine the grid size. Since ^ is defined 
throughout the (q-j, q 2 ) configuration space, the exact grid is 
of infinite area (1 and 3 ujould be infinite). The calculation 
is made feasible by reducing the grid size to a finite range 
and imposing the boundary condition that the u/avef unction is 
zero in all space not encompassed by the mesh. To reduce the 
memory requirement, MW started using L-shaped grid enclosing 
the interaction region, Noiuaday^a rectangular grid is used. 

If the FT method is used to obtain the Laplacian of 
efficiency is maximum if the number of grid points are in 
pou/ers of 2, Typically 64 or 128 grid points are used. The 
FD schemes have no such restrictions on the grid size. As a 
matter of fact, they allouj floating grid size, thus reducing 
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ths memory requirement. For reactive scattering, the Wp at 
time t = 0 must be placed in the entry channel, far out from the 
interaction region (large , small q 2 ). The complex exponen- 
tial term of the Wp [Eq, (52)3 pushes it towards the interaction 
region. The segments of the boundary in both directions which 
close off the asymptotic reactant and product channels consti- 
tute nonphysical barriers which tend to reflect ^ back toward 
the interaction region. Therefore, these boundaries must be 
placed far out (large q-j , q 2 ) so that the reaction is essen- 
tially over before the interference between the waves reflected 
from the inner repulsive wall and the outer boundaries becomes 
significant, 

1,2,3 Calculation of reaction probabilities 

For a specific eigenstate of BC and a distribution of 
R 

, <P > is computed by obtaining the relative volume of 
in the product space P s 

< P^> = f^(t) 5'(ty dq^ dq2 ., (54) 

The products are distinguished from the reactants with the aid 
of a " reasonable ” qj * constant line. For a symmetric exchange 
reaction of the type (R2) the dividing line (of the reactants 
and the products) would pass through the saddle point from the 
origin diagonally* If the saddle point is located asymmetrically 
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in the (q^, q2) space or if there is no saddle point as in the 
case of HeH2 system, the dividing line becomes somewhat arbitrary, 

The initial wavefunction is time evolved for a "suffi- 
ciently long time" and the time evolution is stopped when a plot 
R 

of > as a function of time levels off. Another criterion 
based on the maximisation of entropy principle can be used to 
terminate the time evolution of ^ , The information theoretic 
entropy^^ 

S = -kS Pj^lnP^ = -k / P In P dP 

= -k/ ('?* ^)ln(^* 'P)dq, dq. 

^ ^ (55) 

can be computed at different time intervals and the time evolu- 
tion stopped when a plot of S as a function of time levels off* 
Preliminary results indicate that as time increases, entropy 
increases and then levels off. 

53 

Recently, it has been shown that using the Wp 

f (q^) = expC-ik^q^) sin(Ak( q-]-q-j° ))/( ( q'|-q-)° ) (^Ak) .. (56) 

R 

it is possible to obtain P (E) for a narrow distribution of 
energy. For a WP corresponding to k^^ -k^ and Ak =Ak-| let the 
<P^ be P.|. .for a second Wp with =k2=%^fand Ak = Ak2 = 

Ak^ + e, lct<P^> be P2 . Since Eq. (, 56 > has a rectangular dis- 


tribution in k space. 



over 


R 

ujhere P (E) represents the reaction probability averaged 
the energy values in the range E-AE to E + AE luith 


E = + Ak^ + (58> 

o 

and AE =1^ £(k^ Ak.| + z )/ qq .. (59) 

R 

In principle, AE can be made very small so that P (E) irill 
approach the time-independent result for a fixed energy, 

Leforestier^^ has described how to extract the inelastic 

I p 

(P ^ _) and reactive (P ) state- to- state (m-»n) probabilities 
from the final Wp, For a collisional eigenstate corres- 

ponding to the vibrational state m and to a total energy 


^ C 


m 


8Xp(-ikqT)X^(q2) + S BXp(-ik^q-,)X^( q2) 


,.(k) 


(k) 


^ ) X;/q2) 


.. (60) 


In the above expression, . stands for the dissociative 

^ rn f □! ss 

contribution to that collisional eigenstate and the primed 
symbols refer to the products AB ♦ C. 
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The WP at time t = 0, 

CD 

'^m ^ dk 5 ( k) BXp[-ik( q^-q^° )] .. (6l> 

-CD 

can be projected onto the eigenstates defined by Eq* (60), lead- 
ing to the expression 

\ (qi,^2> = ‘IkJCk) '1^''^ ... (62) 

-03 

The above equation has been obtained by neglecting the 
overlap betiveen the ^' ( k) exp[-ik( q-j-q^° ) 1 function and the 
exp( ikj.^q-|) or axp(ik^iq!|) plane uiaves. This approximation 
becomes exact for high-energy collisions such as CID processes 
because thB$(k) distribution function has essentially no 
components in the negative k plane. 

At time T, the Wp ujould have evolved into 
Vqi,q 2 >T) = J-“clkf(k) exp(-iE T/1f,) .. (63) 

-CD 

If T is sufficiently large, i.e,, if the WP is entirely outside 
the interaction region, one can use the asymptotic evaluation 
of the collisional eigenstates [Eq. (60)1 to reiurite the 

Wp expression as 
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»C|2,T) = -2 ^ dk # ( k) exp(ikj^q^) 

•CX) . ’ 

exp(-iE^ l/iS) 

« , (k^ 

^ dk§(k) exp(ik|!|,qp exp(-iE^ T/ti) 

00 (k) 

dkJ (k) sxp(-iE^ T/f,) ..(64) 

This last expression, after projection onto the and X^, 
eigenfunctions and Fourier transformed yields directly the 

probability amplitudes A^*^^ and 

^ m-*n m-*n * 

(k) -1 

Vn = $(k)] 8xp(iE^ T/1S)(k/k^)* 

x//dq-,dq 2 X^(q 2 ) 8Xp(-ik^q-,) q^ , q 2 ,T) 

(k) -1 

^m^n* “ “[2?rj(k)] Qxp(iE^ tA) (l/'k/A^k|!Ij,) • 

x//dq^ dq^ exp(-ik|l!| , q>|') » T) (65) 


From these equations, one can finally obtain the state- 
to-state probabilities u/hich are defined as the ratio of the 
outgoing flux in the considered channel to the incoming flux. 
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P (k) = h 


rt»*n 




— -l ---^ I/Jdql, dqj X' ( exp(-lk' , q',)? ( q', , qj,T) T 
^k'^, l^(k) ! ^ rf ^ n I m . ^ 

.. ( 66 ) 


1.3 Applications 

1,3,1 Atom-^djatom exchange rsactions 

By far, the predominant application of the TOQM approach 

R 

has been in calculating <P > for a variety of atom»-diatom 
exchange reactions, A summary of them is presented in Taole 1, 
computed the < at different energies for the reaction 
(R2) on the Porter-Karplus^° PES, Agramal and Raff^^ used the 
explicit method to compute rate coefficients in the temperature 
range 1000-3000 K for the collinear reaction (R3) on a London- 
Eyring-Polanyi-Sato (LEPS) PES, The results were shown to be 
in good agreement with the time-independent results of Garrett 
et al,^^ Ualues of < P > computed on the SLTH surface for 
reaction (R2) also have been shown to be in accord with the 
time-independent results, 

Kellerhals et al.^^ computed < P^> for the reaction (R1) 

on a model LEPS surface and studied its dependence Q,n the 

R 

location of the saddle point. It became clear that < P > was 

42 

sensitive to changes in the PES, Stroud et al, investigated 
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TABLE 1, Summary of the atom-diatom exchange reactions studied 
by the TDQm approach, along with a list of numerical 
methods used 


System and 

PES 

Initial conditions 

Laplacian 

Time- 

Ref. 

y 

<^tr> 

evolution 




b\1 




H H 2 


0 

0.14-0.26 

3-point FD 

FD 

27 

model 


1 




H 


0 

0.32-0.51 

3-point FD 

implicit 

29 

PK 







H + H, 


0 

3.00-12.00 

5-point FD 

predictor 

34, 

LEPS 


1 


-corrector 

66 



4 





H + 


0 

2.30-4,10 

7-point FD 

implicit 

35 

LEPS 


2 




H + 


0 

0.40-1.00 

3-point FD 

explicit 

37 

SLTH 


1 





H + 

LEPS 


. 0 

0.15-1.00 

3-point FD 

implicit 

40- 

42, 







62 

H'*‘ +H„; 
DIM 

0 ^* HD« 
D + H^ 

0 

2 


FT ■ 

FD 

44 



3 





T + HD 
LEPS 


0 

0.08-0,21 

3-point FD 

explicit 

53 

T ■*- HD 
SLTH 


0 

0.08-0.87 

3-point FD 

explicit 

56 

F Hot 

F +D_ 

0 

0,50-2.00 

3-point FD 

FD 

61 

Muokerman 






H + H«j 
MOEXP^ 

H +DH 

0 

0.50-16.00 

7-point FD 

explicit 

67, 


1 




6 0 


3 



36 


R 

the influence of v on <P > and its sensitivity to features of 
the PES for the reaction 

He + H 2 ^ HeH'*' + H (R4) 

using a diatomics-in-molecules (DIM) and a slightly different 

1 

spline- fitted ab initio surface. Studies by Zuhrt et al. 
using three different PESs for the reaction 

F + H 2 > FH H (.R51 

R 

also shouJed that <P > was strongly dependent on the topological 

‘ 62 

fFatures of the PES. Sathyamurthy et al, were able to demons- 
trate the utility of cubic splines to interpolation of ^ initio 
PESs by showing that results on a spline- fitted surface were 
identical to those obtained on the original analytic surface. 
Thareja and Sathyamurthy have recently used the TDQM approach 
in testing the quality of an analytic fit to the ^ initia PES 
for the reaction (R2), When different fits were obtained for 
the same original ab initio surface, the fit that led to the 
same <P > as on the SLTH surface was taken to be the best, 

' R 

In addition to yielding P at different levels of resolu- 
tion, the time-evolution method provides valuable insight into 
the detailed nature of the collision dynamics. Typical 3D 
perspective plots of in configuration space at different 
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tiniB intervals give a pictorial representation of the reactants 
transforming into products and the resulting picture can be 
compared directly with a plot of a family of trajectories, 

” St ructur Bs” in these plots reveal the complexity of the dyna- 
mics and are indicative of indirect collisions. Using the FT 

method Kosloff and Kosloff investigated the conditions for 
* 

Cc.mplex formation in H ^2 isotopic analogs for which 

the PES has a deep potential well in it. 

From a knowledge ofW(r,t), quantal flux patterns could 
be drawn, in analogy with the fluid flow in classical hydro- 
dynamics, Plots of the quantal probability flux vector 

=4 I" .. (67) 

29 

were made by MW and quantum whirlpools were identified. These 
have been investigated subsequently by Hirschf elder^^ who des- 
cribed them as quantized vortices around wavef unction nodes, 

1,3,2 Spectroscopy of the transition state 

37 

Recently, we have made use of the TOQM approach to 

predict the wings to the Lyman-a line in H +* H 2 collisions, by 

2 

obtaining the time-averaged If I values and relating them to 
absorption intensity. We have made use of the explicit integra- 
tion scheme^^ to solve the TDSE, Also the same approach has 
been made use of, in predicting the emission spectrum for the 
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system. The results of our study are presented in detail 
in Chapters 2 and 3, respectively. Thus, the approach promises 
to be an effective tool in the emerging area of "spectroscopy 
of the transition state"* 

1,3,3 Collision-induced dissociation processes 

In contrast to the formidable problems faced by the time- 
independent approach in studying CID processes, the time-depen- 
dent approach is applicable to the CID problem u/ith as much 
ease as for the exchange reaction. Ford et al.^^ reported the 
first TDQM results for the reaction of the type 

A + BC > A 8 + C (R6) 

29 

using a model potential and the MW algorithm. The fact that 
the Wp represents an average over a momentum distribution can 
be exploited to extract results for a mide range of energies 
from a single calculation. In this context, the time-dependent 
approach is superior to its counterpart, the time-independent 
approach, Kulander has made use of the TDQM approach in 
computing the dissociation probabilities (P^) for the collinear 
reaction 

H + H 2 ^ H + H + H (RV) 
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on a LEPS surface using the five-point FD scheme to compute the 
kinetic energy term and a fifth order predictor-corrector method 
for the time evolution. 

Subsequently extensive QCT and quantal WP calculations 

66 

have been carried out to compare ujith the earlier reported 
results.^ Leforestier et have studied a variety 

of model CID processes over a ujide range of energies. By 
deconvoluting 4P^> to obtain the energy dependence of for a 
CID process , they have been able to obtain information on the 
threshold behaviour, vibrational enhancement and inhibition and 
their dependence on They have shed further light on the 

dynamical nature of CID processes by plots of 1^1 and 3 at diffe- 
rent time intervals. We have also carried out a study of the CID 
process in the collinear system (R7) on the accurate SLTH PES 
over a luide range of and the results uiill be presented in 

Chapter 4, 

1,3.4 Gas-surface scattering 

Atom-surface scattering plays a vital role in the invest- 
igation of solid surface properties. The rapid advancement in 
the recent past on the experimental side of the atom- surface 
encounters posesa challenge to theoretical models. Time- 
independent quantal calculations have not oecome practicacle 
yet. Classical and semiclassical methods have been used exten- 
sively in simulating the experimontal results. But they are of 
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limited validity. For the first time Agrauial and Raff^® shou/ed 
that the time-dependent approach could be used to study gas- 
aurface scattering by representing an incident atom by a Wp 
and using the explicit integration method^^ for the time evolu- 
tion of the system. Yinnon and Kosloff®^ have been able to 
compute scattering intensities, resonance strengths and other 
attributes using the FT method. Gerber et al. * have used 

the same approach to study the scattering of a He atom from a 

72 

Cu surface luith isolated Ar impurities, Kosloff and Cerjan 
have studied the desorption and scattering of a He atom from 
W and Pt surfaces using the same method, 

73 

Recently, Moujrey and Kouri have shomn that by solving 
the TDSE, molecule-surface scattering processes can also be 
studied. They have investigated inelastic scattering of H 2 from 
a flat, rigid surface using a hybrid scheme combining a coupled- 
channel expansion u/ith a Wp formalism. 

1 .4 Highlights 

The TDQM approach scores over the time-independent approach 
in certain respects, in that there is no need for the construc- 
tion of channel hamiltonians and matching of the ivavef unctions. 
There is no increase in the complexity of the problem or the 
computational time reguirement udth an increase in the total 
energy of the system. As a matter of fact, there is a reduccion 



41 


in the computer time uuith an increase in energy-analogous to 

the classical approach. In contrast to the time-independent 

R 

approach in ujhich P has to be computed for different and 

\i and averaged over the appropriate distributions in order to 

calculate state-selected and overall rate constants, the TDQn 

method could be made to yield them directly by a suitable choice 

of initial ^ . Also, the deconvolution methods have made it 

R D 

possible to obtain in a single calculation, P (or P ) over a 
uiide range of 
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WINGS TO THE LYMAN-a LINE IN H+H 2 COLLISIONS 

7 4-80 

It has been recognised in recent years that it is 

possible to obtain the spectrum of the transition state (me 

use this term in a general sense to refer to all configurations 

intermediate betmeen reactants and products) for elementary 

7 5 76 

chemical reactions. Arromsmith et al, * reported the 
emission spectrum of the TS for the reaction 

Ma^ + F > NaF + Na* 

by monitoring the "mings" to the sodium D-line emission,. Polanyi 

7 7 7 R 

and Wolf * interpreted the wings as arising from the emission 
from Na* while the separating NaF is still in the v/icinity. In 
other words, they concluded that [FNaNa] was the emitter. The 
emission signals were observed upto~ 400 S away from the sodium 
D-line, in spectral regions where no transitions are known for 
either reactants or products. Pressures were orders-of-magnit uo- 
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too low for the signal to arise from far wing collisional 
broadening and hence the observed emission must have been due 
to the TS, 

Brooks and coworker s"^^ have reported evidence for absorp- 
tion from the intermediate configurations in the reaction 

K + HgBr 2 ^ KBr + HgBr (R9) 

Although K and HgBr 2 react on nearly every single collision to 
form KBr and HgBr (X 2 ) , the exoergicity is insufficient to 
yield the electronically excited HgBr*. However, if the system 
absorbs a photon with x 590 nm [no single photon absorptions 
are known for either reactants or products at this wavelength] 
the electronically excited HgBr* (B ) state, becomes access! 
ble. And this state is easily detectable from the fluorescence 
(§ — 5^), Thus, it was shown that (R9) is a laser-driven 
chemical reaction in which the excitation of the reaction 
complex opens a new channel 

K ♦ HgBr 2 + hy (= 595 nm) KBr + HgBr* (i (RIO) 

Fluore- 
scence 
500 nm 

HgBr (X '^Z*) 

80 

Recently, Brooks and coworkers have also reported evidence 
for light absorption by the reaction complex in a similar 
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reaction 


K + NaCl + hv (v^ 700 nm) — » KCl + Na* (3 ^P) 

Emission 
V 5 89 nm 

Na (3 ^S) 

81 8 ? 

The studies * on the reaction 


(R11) 


Xe + CI 2 + hv XeCl*" + Cl (R12> 

also indicate the reaction complex absorbing light, 

03 

Spectra of such TSs can provide information on the 

collision dynamics, complementary to vuhat can be learned from 

84 

state- to- state studies. 

It is of fundamental interest to predict and observe the 

u/ings to the Lyman-a line in H + H 2 collisions. Polanyi and 
85 86 

couiorkers * have computed the absorption spectrum for the 
TS for the collinear reaction (R2) by determining the probabi- 
lity density distribution for the TS on the chemically accurate 
SLTH surface using the QCT technique. They assumed vertical 
transitions from the ground state to one of the loUJ lying excited 
states, ^TT in geometry , that correlates mith H*(2p) ♦ H 2 . 

For a quick reference, a correlation diagram of a feu/ louj lying 
states of H 3 is presented in Fig. 5. They used model PESs based 
on the auailable ab initio data for the upper PES. They studied 
the effect of increasing relative translational energy of 




88 

correlation energy diagram for the first few lowest 
potential- energy surfaces of together with the 

ground state surface as a function of nuclear geometry 
The lowest excited-state ^d the ground-state surfaces 
meet in a conical intersection at D 3 . geometry 


Pig. 5. 
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reactants and reagent vibrational excitation on the features 
of the absorption spectrum. As the collision energy was 
increased, the intensity of the spectrum diminished and there 
was a considerable red-shift for the wings to the Lyman-ix line. 
Subsequently, they extended®^ their study to collisions in 3D 

Qq 

using a diatomics-in-molecules (DIM) surface for the excited 
state and concluded that the qualitative features of the 
spectrum remained while going from collinear geometry to 3D, 

90 

Recently, Engel et al. have carried out time-independent 
QM calculations and predicted the absorption spectrum of the TS 
for the collinear reaction (R2) for two different vibrational 
states of the H 2 molecule, at collision energies below the 
reaction threshold (0,22 bU) using the same ground and upper- 
state surfaces as Polanyi arid coworkers did. They computed 
the continuum- bound absorption cross section using the relation 

a (v) I 2 .. (68) 

” tr 

where denotes the excited state wavefunction describing the 

n— th bound state and , the continuum wavefunction at 

tr 

They treated the transition moment ^J to be a constant (=1) as 
it was done in earlier studies, * The bound . states, covering 
75 per cent of the potential well depth on the upper PES were 
calculated by expanding the wavefunction in terms of an ortho^^ 
normal set of vibrational eigenfunctions. The continuum 
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ujavefunction u/as calculated using three vibrational expansion 
functions. Their generated spectra agree well with the spectra 
computed via the QCT method with respect to the basic trends 
and features. The overall appearance of the absorption spectrum 
has been interpreted in terms of the amplitude pattern of the 
continuum wavefunction and the radial extent of the excited 

state bound level wavefunction. Subsequently they have extended 

9 1 

their study to reactive H + H 2 collisions. 

92 

Recently a model study has been carried out in our 

laboratory in which the absorption spectrum for the TS in the 

collinear reaction (R2) has been predicted, using the time- 

independent SchrQdinger equation for a model consisting of a 

ID rectangular potential barrier and a well, at three different 

collision energies. Results have been shown to be in qualita- 

9 0 

tivB agreement with those obtained from a full QM study, 

37 

We have also computed the absorption spectrum for the 
TS for the collinear reaction (R2) , but using a TDQM approach. 

We have used the SLTH^ PES for the ground state and a DIM 
surface (described in ref, 88) for the excited state, that 
correlates with H*(2p) + H 2 . The contours of the ground and 
the excited state PESs are displayed in Fig. 6 and Fig. 7, 
respectively. The TS spectrum has been predicted for four 
dif ferent <E. >. [two of the energies lie below the reaction 

threshold and the other two well above the threshold] for two 
different vibrational states (v=0 and 1) of the H 2 molBculB, 
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TS spectra have also been computed for thermal distributions 
of H + H 2 , at 300 and 1000 K. 

2.1 Computational detaj 1 s 

The tujo dimensional TDSE has been solved by using the 

53 37 

explicit integration method. In our earlier work, we had 
used a 3-.point FD scheme to compute the kinetic energy term in 
Eq » ( 59 ). Subsequently, 5-point. and 7-point FD schemes and 

also the FFT algorithm have been used to compute the same. The 
results are summarised below, A (59 xA5) rectangular grid in 
the (q..|,q 2 ) space with 1 ,7 < q..| < 8,08 a.u. and 0,4 ^q 2 < 5,24 a.u. 
has been used. The numerical values for units of distance, 
time, mass and energy used in our calculations are given in 
Table 2, The initial WP was located at q-j° = 5,5 a.u, and had 
a width of 6 = 0,25 a.u. The coordinate mesh spacing ifisq-| = 

^q 2 *-0,11 a.u, and the time step At = 0,01 were chosen such that 
they meet the stability criteria for the explicit integration 
method as discussed in Chapter 1, The initial wavefunction 
f(t *0) was set up as in Eqs, (46) and (52) for different 
of interest. 

The time evolution of the system described by TDSE was 

R 

terminated at a time when the plot of <P > as a function of 
time levelled off as shown in Fig. 8, We have checked the 
accuracy of our calculation by checking the normalisation 
Condition for 1 throughout the time evolution, of the system. 



/Vo,. 
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TABLE 2. Molecular Units 


Quantity 


Mass 

Distance 

Energy 


1 atomic mass unit = 1 
1 atomic unit = 1 a.u 
1 electron volt = 1 eW 


1 t.u 


a.m.u. = 1.6604345 x 10 
. = 0.529167 X 10“® cm 
= 1.60210 X 10“^^ erg 
. = 0.5 3 87 1 47 X 10“^^ 


-24 


g 


Time 


1 time unit 


sec 


REACTION PROBABILITY 
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Fig. 8. Reaction probability as a function of time 
with each time step = 0.53875 x 10 s, f' 
a spBcifiBd initial condition 
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Ths accuracy of ths method has also been checked by back-evolu- 
tion technique. That is, the WP is first time evolved from 
time to to t^ + nAt, Then, by replacing At in Eq. (21) by -At, 
the u/avefunction is "back evolved" n steps to time t^. Using 
this technique, uie luere able to reproduce the initial wave— 
function accurately* If the computed <P > is a converged one, 
then it should be independent of the number of grid points used 
and also on the boundary conditions* We have verified that the 
above statement is true, by computing <P > uiith a square compu-* 
tational grid (64 x64) space (q^,q 2 ) ujith T^T^q^ <11^15 a*u* 

and 0, 4 S. 9 • 85 a*u* A more stringent test on the accuracy 

R 

of our calculations mould be to compare our computed <P > to an 

R 

average of the time-independent reaction probability P (E) over 
the distribution used for the WP, in the folloming may. 


o /P(E) P^(E) dE 

<p > = — 

/P(E) dE 


.. (69) 


mhsre the meighting function P(E) is given by 

p(E) = l#(k)l (7 0) 

Substituting for#(k) and using the relations 
k = (2mE)^^^/fi 

dk “ -(2m)^^^M]E:~^'^^ clE .. (7 1) 
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We get 

<E>= (26^/n)'’/2(m/2)^'^^(l/1S)exp[-26^(2m/4^^)(E^'^^-Ey^)^] 

o 

xE-‘’/^.E.dE ..(72) 

Thus, 

P(E) = (26V7t)''/2(n,/2)''/2(i/if,). 

axp[-262(2m/lS^')(E^/2^£V2^2j^-l/2 

R 

Substituting Eq, (73) and the values of P (E) reported in ref. 60 

R 

into Eq, (69) and performing the integration, me get the <P > 
that is to be compared with our computed results. The results 
of these comparative studies are presented in Table 3, Since 
me have computed <P > for different WPs that correspond to 
Maxmell-Boltzmann like distribution of relative speeds at tempe- 
ratures of 300 and 1000 K as mell, the corresponding rate 
coefficients (k) can also be calculated using the relationship 

k = (kgT/2^ii)^^^ <P^> ..(7^> 

(mhere kg = Boltzmann constant) and compared mith the results 

^ ^ 60 

reported in the literature. 

We have tested the accuracy of our results obtained by 
using the 3-point FD scheme to compute the Laplacian of 'F by 
repeating the calculations using 5-polnt and 7-point FD eohemes 
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TABLE 3. Comparison of reaction probabilitiss and rata coeffi- 
cients for the collinear reaction H +“ H 2 H 2 + H as 
computed from the TDQM approach and its time-indepen- 
dent analog 


V 

Temp 

(K) 

^^tr ^ 

nR 

<P > 

<p'’> 

k ( 10 ‘^cm 

/ sec ) 


(eV) 

(a) 

(b) 

(a) 

(b) 

0 


0.41 

0.561 

0,637 

- 

- 

0 

- 

1.00 

0.497 

0.490 

- 

- 

0 

1000 

- 

- 

- 

0.61 

0.67 

1 

1000 

- 

- 

mm 

2.71 

2.98 

(a) 

present 

calculations 






(b) ref. 60 
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and ths fT mathod based on the FFT algorithm.^^ In the 

latter cases^ norni for ^ was not conserved when we used the 

value of Aq = 0,11 a.u, as was used in the 3-.point FD scheme. 

Therefore, we had to use a larger value of Aq = 0,15 a.u. In 

the case of the higher-point FD schemes and the FT method, use 

of a small value for mesh spacing makes the round-off and the 

truncation errors dominate and hence spurious results are 

prottiuced. So the argument that use of a smaller mesh spacing 

always makes the results to^improve, becomes invalid in this 

context. In fact, only above a critical value of the mesh 

spacing these formulae for second derivatives become operative 

R 

and effective as well. We got almost identical results of<P > 
in all oases, as it is illustrated in Fig, 9. 

For a specified set of initial conditions, v =0 of the 

H 2 molecule and = 0,41 eU with a 64 x64 grid, the mesh 

Space being 0,15 a.u,, a run involving the 3-point FO scheme 

took 12 min of CPU time on the DEC-IO 9 O computer. The corres- 

32 

ponding run involving the FT method took 180 minutes. As has 
aiready been mentioned in Chapter 1, the number of operations 
per grid point is larger for the FT method and thus it is 
15 times slower than the FD approach. 

We solved the TOSE for H + H 2 motion on the SLTH surface 
as described above and obtained lP(t> at different time inter- 
vals. Then we adopted the following strategy to compute the 
TS spectrum. The value of l?*(t) W (t) dq^dq 2 over each area 
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TIME 


Fig. 9. Reaction probability as a function of time, 
using different schemes to compute the 
Laplacian, for a test case 
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element is taken to be the probability density for the TS in 
that area element. The time-averaged value of W*(t) 'I'(t)dq^dq 2 

over the entire time evolution period for each area 
element represents the intensity of the transition at an energy 

equal to the difference between the ground and the excited- 

state surfaces in that area element. We have used the DIM 

Surface described in ref, 88 for the upper surface. We have 

assumed constant transition moments as was done by the earlier 
86 88 90 

workers, * ’ We have used a histogram width of 0.2 eU to 

compute the spectrum. 

2,2 Absorption Spectrum for the TS 

The resulting absorption spectra for \^^ for rang- 

ing from 0,10 to 1,00 aU for v»0 and 1 are displayed in Fig. 10, 
It is clear that for high energy collisions there are substan- 
tial wings to the Lyman-a line at the lower energy region. Two 
prominent peaks have been identified (at 72,600 and 51,600 cm ) 
as a and b. With an increase in the collision energy, the 
intensity of the peak close to the line decreases and that 

of the peak in the lower frequency region increases. Similar 

86,88,90 

trends have bsan report ed by earlier uiorkers# 

These results can be understood as follows. We show the 
comulative value of the probability density function (PDF) 
over the entire grid, superimposed on the difference potential 
contours (AU ^ difference in the potential energy 


Intensity in arbitrary units 


Fig. 



10. Absorption spectra for under different rnrtxal 

conditions. Average collision energies are indicated 
against each spectrum. represents the Lyman-a 

line. 
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between the upper and the lower surfaces) in Fig. 11 and 12 for 
d».ff8rent initial conditions. Hence each point in the confi- 
guration Space can be read as a spectral frequency v, corres- 
ponding to the energy difference AU. The greater the PDF at 
any particular AV , the greater the absorption intensity at the 
corresponding v. Thus, the peaks in the absorption spectra can 
b© identified with the accumulation of the PDF at a particular 
range of AV values. At low collision energies (0,10 and 0,25 e\l) 
do not SCO substantial wings to the L line. At higher 
Ofttllision onorgios (0,41 and 1.00 sV) two peaks in the lower 

BfiBrgy (frequency) regions have been identified as a 

1 *1 

(“?1 ,600 cm* ) and b (51,600 cm" ), These peaks reflect the 
build up of the PDF near the inner repulsive wall in the PES, 
Examination of Fig, 11 shows that at low PDF is con- 

f ined to the entry channel indicating that the Wp is not able 
to cross the barrier and hence retrieves back to the reactant 
valley. As <E^j,>is increased, PDF is smeared all over the 
c onfiguration space and there is a suostantial build up of PDF in 
particular regions of space. The peak labeled a in Fig. 10 
f-or identified with the probability 

density contour (PDC) marked 3 and 5 running parallel to the 
aixis in the region between 8 and 9 eV, The peak labelled 
b is duo to tho build up of PDF marked 3 near the saddle point. 
Pis the collision energy is increased to 1.00 eU , the PDC marked 
Sis found to bo missing in the AU region between 8 and 9 oM 
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and hence the intensity of the peak " a “ falls off. But there 
is a build up of the PDF narked 4 in the AV region betiueen 5 ^ 6 bM 
explaining the increase in the intensity of the peak "b" with 
increase in . 

Within the resolution of our spectra, we do not see any 
significant effect of vibrational effect of H 2 on the spectrum 
of the TS, 

Transition state spectra, computed for thermal distribu- 
tions of H H 2 at 300 and 1000 K are displayed in Fig, 13, In 
contrast to the spectra for fixed <E. > (in Fig, 11-12), the 

spectra in Tig, 13 do not show any significant "far” wing as the 
<£tjj.^ corresponding to T *300 and 1000 K are only 0.03 and 
0,09 ey. The PDF plots in Fig, 14 show that the Wp is confined 
to the entry valley and does not reach upto the inner repulsive 
wall because of the low >, This would explain why the 

"wings" ware not discovered long ago, although atomic spectra 
have been recorded for a long time. We must add that the 
shoulder to the Lyman-a line becomes marked with increase in 
T from 300 to 1000 K in keeping with the observations made for 
the "fixed" results. 

If the complementary information regarding the emission 
spectrum becomes available, the analysis of the TS spectrum 
would become easier. We have computed the emission spectrum of the 
TS for the same reaction (R2) by carrying out the dynamics on 
the PES for the excited state of that correlates with 



Intensity in arbitrary units 


Y = 1 


f 

I 

I 



g. 13. Absorijtion spectra for nf for a thermal distri- 
bution of H, H 2 (v =0 and l) at 300 and 1000 K. 
The dashed line indicates the spectrum obtained 
in the absence of any interaction between H and 

(at t = 0) 



Fig. 14. Same as Fig. 11 but under dirferen 
as Indicated in each panel 
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H*(?p) "^^2* Cistails of such a calculation and the results 
are reported in Chapter 3, The effect of including non- 
collinear geometries, a variable transition moment and tran- 
sitions to other lou)-lying upper electronic states and its 
implications to an experimental observation of the wings to 
the Lyman-cr line are currently under investigation. 



CHAPTEfi 3 


EMISSION SPECTRUM OF THE TRANSITION STATE wf* 

By recording the emission spectrum of the transition state 
for a chemical reaction, one can get a oiealth of information 
about the dynamics of the system. From a knowledge of the total 
integrated intensity, the breadth and the spectral distribution 
of the wings, one can learn a great deal regarding the lifetime 
of the TS, preferred turning points on the PES, bending in the 
TS and attractive or repulsive character of the PES, That such 
a link exists between these factors have been demonstrated by 
preliminary theoretical investigations, * The TS in a reac- 
tion of the type (R1) reveals its existence and the time it 
spends at successive configurations by emitting radiation of 
corresponding frequency and intensity as illustrated schematic- 
ally in Fig, 15, 


Emission spectra of dissociating molecules appear to bo 
particularly promising for obtaining detailed information on 
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PRODUCTS SEPARATE 


Fig. 15. Schematic diagram of emission from the TS 
in a typical chemical reaction.- Emission 
from ABC during separation of products 
results in "wings" to the atomic transition 
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the TS for the corresponding bimolecular chemical reaction in 
the reverse direction. 

"I O 

It has been sho^n that by recording the photoemission 
spBctrum of (dissociating) polyatomic molecules^ the geometry 
of th© upper state PES near its saddle point could be determined 
and hence the possibility of spectroscopic probing of dynamics 
has been demonstrated. By his SC approach, Halier'^ ’ has 
shown that such problems are amenable to theoretical treatment 
at present. In this approach, the ground-state initial wave- 
function (at t * 0) is transferred to the electronically excited 
stoto and then the time evolution of the Wp is governed by the 
nucloar-motion Hamiltonian whose potential is that of the excited 
state, A simple Fourier transform in the time domain brings us 
to the frequency domain, Heller has shown that it is usually 
unnecessary to follow the complete evolution of the WP on the 
upper state and that short-time behaviour can be described 
adoquatoly by slmplo SC dynamics. He also derived simple 
formulae to make it possible to infer the characteristics of 
the upper-state PCS directly from the emission spectrum. 

Imre and coworkera^^ have recorded the spectrum of 0* 
molecule. Conventional spectroscopic techniques could not be 
applied to study the excited states of 0* owing to their extre- 
mely short lifetimes. Some of the elBctronically excited states 
are not bound, thus leading to rapid dissociation. From 
an experimental knowledge of relative intensities of the 
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fundamantal symmetric stretch (j;^) and the first overtone of 

the asymmetric stretch (vj), they have been able to extract 

quantitative information regarding the excited state under 

study. By inserting their data in the relevant equation derived 
12 a 

by Heller, they obtained a negative value for the second 
dorivativo of the potential with respect to the Qj normal co- 
ordinate, indicating that the upper PES has a maximum at the 
ground-state equilibrium geometry. This result has been vindi- 
cated subsequently by initio computations.^^ 

In an attempt to complement our studies (Chapter 2) on 
the absorption spectrum of the TS for the reaction (R2) , 
we have carried out a TDQM study to predict the emission 
spectrum of the TS for the same reaction. We have studied the 
dynamics on the PES for one of the low-lying excited states of 
Hj that correlates with H*(2p) H 2 and considered the emission 
to the ground state. The same pair of PESs that were used in 
predicting the absorption spectrum in Chapter 2 have been used 
in the omission study, 

’fc 

The potential energy contour of the Hj PES has been dis- 
played in Fig, 7, The initial wauefunction 'P(t=0) has been 
set up in the seme manner as was described in Chapter 2. A 
64 x64 grid in the (qT,q 2 ) space with 1.0 < q^ <7.3 a.u. and 
0,4s,q2i.6.7 a.u. has been used. Attempts to use 3-point and 
5-point FD schemes with Aq=0,1 a.u. and At = 0.01 (the values 

that have been used in the absorption spectrum caloulations) 
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failed. The norm in f mas not conserx/ed. Since the PES has a 
well of about 4 eV [see Fig. ?], in a small range of the inter- 
action region, we faced problems regarding the numerical stabi- 
lity of the solution. Therefore we had to choose a prohibi- 
tively small step-size in time, At = 0.QQ1. The 7-point FD 
scheme with a mesh spacing of Aq = 0.1 a.u. had been used to 
compute the Laplacian of I? . We found that the norm in ^ was 
conserved throughout the collision for the choice of these 
values of At and Aq. We are not surprised that for the choice 
of a small value of Aq, the 7-point FD scheme was stable (in 
contrast to the report we have made in Chapter 2), The diffe- 
rence, however, was due to the choice of a small value for At. 

An important factor in integrating partial differential equa- 
tions [£q,(1)] is the ratio (Aq/At) and not Aq or At alone. 
Accordingly the ratio of Aq and At was adjusted and the devia- 
tion of the norm of f from unity was closely followed, to ensure 
the accuracy of the calculation. 

We have considered only collisions with the H 2 molecule 
in its ground vibrational state and 

used the explicit method (as described in Chapter 2) to time 
evolve * and the results are shown in the form of contours of 
at different time intervals in Fig, 16. It is interesting 
to note that the Wp tends to spend considerable amount of time 
in the region of the configuration space where the well is 
located [panel (b) of Fig. 16], We also obtain a <p''> of 0.4 



Qj/bohr 


72 



Pig, 16, The probability density of the system at different 
time intervals. The dashed line indicates the zero 
kcal/mole PE contour. The hatched region in panel (a) 
indicates the probability density at time t=0. The 
probauility density at different time intervals 
(indicated against each panel) has been shown by 
solid lines 
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in reasonajlG agroement mith the quasiclassical trajectory 

94.. 

result obtained under identical conditions. At the end of 

the t i fTi 0 evolution of the system^ the probability density 

becomes smoarod all over the configuration space [panel (d) in 

Fig. ^6"^. The same feature has also been observed in the QGT 
94 

results. 

Using the procedure described in Chapter 2, ujb have 
generated the emission spectrum of the TS and the result is 
displayed in Fig. 17, Three prominent wings labelled "a" 

(77,500 cm”S, ” b*' (75,000 cm"*') and "c" (64,000 cm"*”) to the 

Lyman-a line have bean predicted. In addition, in the far wing 

-1 

region n hump *'d'* around 52,500 cm is also noticeable. Since 
the PDF distribution superimposed on the A\i contours, shown in 
Fig. 18 is so complicated and since we have the results only 
for a single collision energy, we have not been able to clearly 
identify the wings with the build up of PDF in the corresponding 
region of tho configuration space as wa could do for the absorp- 
tion sp* Ctrurn. A moro detailed study to understand the structure 
of tho spoctrum is in progress. 









CHAPTER 4 


COLLISION-INDUCED DISSOCIATION IN C.OLLINEAR 
H + H 2 COLLISIONS 


Collision-induced dissociation processes of the type 

A + BC — > A + B + C (R6) 

play a vital role in the kinetics of high-temperature systems 
such as shock tubes and flames. They also have a special 
relevance in connection uJith unimolecular decay and three- 
body recombination processes. CID is also one of the many 
possible outcomes in high-energy crossed molecular beam reac- 
tions, electron attachment and charge transfer processes. 

Considerable amount of experimental results on systems 

95^9 8 

exhibiting CID has been available for some time. Because 

of the difficulties associated in computing CID cross sections 
this process has received less attention from theorists irhen 
compared to elastic, inelastic and exchange reaction processes 



77 


The available theoretical results on CID processes are based 

on phase-space theory, optical model, qQ.j-l02-.107 
1 0 110 

SC» Results of a feu; time-independent QM calculations 

restricted to collinear geometries have also been report- 
ed, * “ The major difficulty in the latter approach is 

the presence of the dissociative continuum preventing a 
straightforu/ard extension of the methods developed for study- 
ing inelastic and exchange reactive scattering processes. A 
comprehensive revieui of the close-coupling and QCT methods for 
CID calculations is given elseuihere, ^ ^ ^ 

As has been pointed out in Chapter 1 of this thesis, the 
TDQM approach can be used to study CID processes with as much 
ease as for the exchange reaction, Very recently, the TDQM 
method has been used to study CID processes in several 
systems.^^*^^ Although they have been restricted to 
collinear collisions, they have provided an insight into the 
dynamics of CIO processes and have helped in assessing the 
validity of approximate methods. 

In real systems, at high energies, the CID process may 
compete with an exchange process. An accurate quantal method 
should be able to incorporate such a possibility. The develop 
ment of such accurate methods have been hindered by the diffi- 
culty of representing exchange product bound states in terms 

111114 

of the reagent bound and continuum states, * As a result 

mostly systems in which exchange channel is absent have been 
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considered in previous studies. Only recently, 
studies of systems in which exchange reaction may also occur, 
have been carried out.^ 36a yj^g TDQM msthod was first applied 
to CID processes by Ford et al,^^ using a model potential. In 
the TDQM approach the problem associated with a dissociative 
continuum is circumvented by a completely numerical descrip-. 
tion of the wavefunction and its time evolution. The disadvan- 
tage of the TOQM method is that it requires the solution of the 
TDSE for evolution of the wavefunction over a long period of 
time corresponding to the entire collision event. With the 
advent of fast computing machines, this is not a serious 
problem anymore. 

34a 

Kulander reported the results of a collinear CID 
calculation for the first time, employing a realistic PFS that 
corresponded to (H , H 2 ) system. He solved the TDSE for the 
chosen Wps and projected the time-evolved wavefunction onto 
the asymptotic eigenstates of the system to extract state-to- 
state transition probabilities. He used a 5-point FD scheme 
to compute the Laplacian of ^ and a fifth order predictor- 
corrector method to time evolve the system. No exchange 
reaction was observed around the energetic threshold ) for 

dissociation. The dynamical threshold ) for dissocia- 

tion was substantially larger than the E®j^ for the ground (v-D) 
and first excited (v=1) state of the molecule. However, when 
the H 2 molecule was excited to v = 4 state, substantial 
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dissociation occurred. Gray et observed no dissocia- 

tion for V =0 and 1 states of the H 2 molecule, when the calcu- 
lations were carried out using the QCT method for the same 
system under identical conditions. Subsequently the quantal 
calculations were repeated^^ to compare quantitatively the 
quantal results with the QCT results^^ and the same trend was 
observed. Quantal was only 0.1 for v=0 and 1 states of the 
molecule at as high as 12 eU, QCT results snowed no 

dissociation at all upto the <E^^> considered, 

67 

Leforestier et al, have carried out a TDQM study in 

1 10 

order to assess the validity of the SC approach to CIO pro- 
cesses. Their system consisted of a Morse oscillator inter- 
acting through an exponential repulsive potential (MOEXp) and 
hence did not allow the possibility of an exchange reaction. 

They used a 7-point FD scheme to compute the Laplacian of ? and 

31 

the explicit method for the time evolution of the system. They 

found the E^^ to be twice the More recently, they have 

68 D 

extended their study to compute P with the molecule under 

attack in different initial vibrational states. In this study 

also, they observed the E^^ to be twice the E®^ . They attri- 

1 07 

buted it to the collinearity of the model, 3D QCT results 

suggest that the CID process takes place mainly by non-colli- 

near collisions. The good agreement between quantal and the 

67 6 8 

SC dissociation probabilities showed * that the latter 
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had carried out a TOQM study of the CIO process in (H,HD) 
system where inelastic and exchange scattering also occurred 
simultaneously. He used a modified LEPS surface; a 7-point 
FD scheme to compute the Laplacian of f and the MW algorithm^^ 
for the time evolution of the system. Both the and the 

coincided with one another. Earlier wor|<ers^^^^° had also 
obtained a similar result. This effect is attributed to the 
fact that energy transfer is more efficient in reactive colli- 
sions than in inelastic ones. 

116 

Recent 3D QCT calculations for three different systems 
(HtHj) , (Hb,H 2 * ) and (He ,H 2 ) on accurate _ab initio PESs 
indicate the near coincidence of the E^^ for dissociation with 
that of the E®|^ for two of the three systems. While (H , H 2 ) 
and (He,H 2 '*^) are potentially reactive systems, (He, H 2 ) system 
is not. Hence for the latter system , The basic 

difference between the results for the different systems have 
been attributed to the difference in the shape of the PE3 
involved, A more comprehensive study is needed to undenstand 
the factors governing CID processes. 

We have carried out a TDQM study of the CID process in 
collinear (H, H 2 ) system over a wide range of energy. We have 
used the chemically accurate SLTH surface for the system. The 
molecule has been represented by a Morse oscillator and the 
classical dissociation energy has a value of ^.75 eU, We have 
used a 7-point FD scheme to compute the Laplacian and the 
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Explicit method to time-evolue the system. Because of the 

high collision energies involved, lue had to use a small step- 

size in time, iit =0.001. A 64 x64 computational grid in 

space ujith 1.0<q^<7,3 a.u. and 0.4 <q 2 <6.7 a. u. 

and Aq-| = Aq 2 * 0,1 a.u, has been used. At energies well 

celou) the dissociation threshold me observed substantial 

exchange reaction probability which became zero around 

When, the mols’cul© was chosen in its ground vibrational state 

(v=0), no dissociation was found to occur even at <E^^> 

energies as high as 10.0 eV, But when the molecule had been 

chosen initially in its first excited state (v=1), at the 

corresponding the dissociation probability rose to 

around ten per cent. Thus, the has been found to be twice 

the E®|^ for the first vibrational state, in accord with earlier 
3D6 

reported results for the same system but on an LEPS surface 
using the QCT method. Our results thus agree qualitatively 
with other quantal result 

In confirmity with earlier explanations, we also find 
that the difference in the dynamical and energetic thresholds 
arises mainly from the constraint that the three nuclei lie 
along a straight line. Also, this feature has been attributed 
to the fact that the exchange reaction probability is identi- 
cally zero at E®J^, We have also observed this behaviour as is 
illustrated in Fig. 19. Fig. 20 shows the probaoility densitie^ 
at different time intervals during the evolution, for v=0 stat£| 
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19. Exchange and dissociation probabilities as a 

function of collision energy. The arrows indicate 
the energetic threshold for dissociation. Note 
that P = 0.0 for v = 0 of H 2 molecule. The 
hatched portion indicates an upper bound for 
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Fig. 20. Probability densitiea in the configuration space at 
different time intervals for v =0, <E^j^>= 10.0 eV . 

The dotted lines are PE contours with values indicated 
in eV against them. The solid lines are the contours 
of X indicates the saddle point. 
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of the molecule at <E^^> = 10,0 eU, It is evident that there 
is a large amount of elastic scattering* A close inspection 
of fig. 20 reveals that the Wp is reflected off the repulsive 
wall of the PES, In Fig, 21, the time evolution of the system 
with the molecule in its first vibrationally excited state has 
been depicted. The node in the initial ujavef unction is 
retained even at the end of the collision. Moreover, addi- 
tional nodal structures become apparent, suggesting a propen- 
sity rule for vibrational excitation, dv = 2, This observa- 
tion can be verified if we project the final oiavef unction onto 
the asymptotic eigenstates of the reagents and compute the 
stat e-to-stat e transition probabilities, A more detailed study 
is in progress. 



Fig. 21. Same as Fig. 20 , for v = 1, < > = 10.0 eV. The 

dashed line indicates the nodal structure in the 


wavef unction 







CHAPTER 5 


SUMMARY AMD CONCLUSION 


We have revieiued the different computational methods 
that are available till to-date to carry out a study of 
reactive scattering processes using the TDQM approach and 
their applications to problems of chemical interest. We have 
carried out a comparative study of the different algorithms to 
compute the Laplacian of the uiavef unction. Even though the FT 
method makes use of the accurate nonlocal approximation, it 
has been shown that it is computationally fifteen times slower 
than the FD schemes, for the same level of accuracy. 

We have predicted the absorption spectrum of the TS for 
the collinear reaction (R2) by solving the TDSE for nuclear 
motion on the chemically accurate SLTH PES, We have made use 
oT the explicit method for the time evolution of the system. 

For the upper surface, we have used a DIM PES , that 
corresponds to H*(2p)+H2. Absorption from the lower (SLTH) PES 
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to this surface been predicted, assuming constant transition 
moment. The u/ings to the Lyman-a line have been interpreted 
in terms of the build up of the PDF near the inner repulsive 
u^all of the PES. The effect of varying relative translational 
energy of the reactants and the vibrational state of the mole- 
cule have been investigated. As the collision energy is 
increased, the intensity of the luing close to the Lyman-a line 
falls off and the intensity of the far wing increases, \/ibra- 
tional excitation of the molecule does not affect the features 

of the spectrum. Our results are in agreement u/ith those 

9D 

obtained by using the time-independent Qfl approach as uell as the 
86 8 8 

QCT ' approach. We have computed the TS spectra for thermal 
distributions of H + H 2 » at 300 and 1000 K, We did not observe 
any ivings to the Lyman-a line under thermal conditions as 
that corresponds to these temperatures is too louj to give rise 
to observable luings. 

We have carried out a dynamical study on the PES which 
has a potential luell of eW using the TDQM approach. A 7-point 
FD scheme to evaluate the Laplacian and the explicit method 
for time evolution of the system have been used. By generating 
the PDF on this PES, we have predicted the emission spectrum 
of the TS and it is shown to be complementary to the absorption 
spectrum results reported in Chapter 2, 

We have investigated the CID process in the collinear 
H +H 2 collisions (R6) in detail, on the chemiGally accurate 
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SLTH surface, using the TDQM approach. We have computed the 
exchange and dissociation probabilities over a u/ide range of 

(0.5-10.0 eU) . No dissociation has been found to occur 
uiith the molecule in its ground vibrational state (v=0) even at 
the highest ^ ^tr) studied. For v=1 state of the molecule, 
has been found to be greater than tujice the This feature, 

we have attributed (in conformity with the results of earlier 
Workers) to the fact that the exchange reaction probability is 
zero around for dissociation. Both the dynamical and ener- 

getic thresholds have been reported in the literature to co- 
incide with each other for systems wherein the dissociative 
process is accompanied by the reactive (exchange) process. This 
feature had been explained in terms of the efficient transfer 
of energy in reactive in reactive rather than inelastic collisions 

Our results are in agreement with other results computed by 
using TDQM^^»^^ and QCT^°^ methods. Recent 3D QCT^^^ results (on : 
the same PES) for the same system indicate that both the thres- 
holds and E^^) agree well with each other. Since experi- ' 

96 1 07 

mental results^ and other 30 QCT calculations on similar 
systems have also shown that both these thresholds usually match,! 
we have rationalised the difference between the two thresholds 
we have observed to be an artefact of the collinear model. In 
addition, careful inspection of the PDF plots of the system with 
the molecule in its v = 1 state revealed new nodal structures, | 

in addition to the node present in the initial wavefunctlon, | 

This fact suggests that this collisional process may be gauerned b 
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a propensity rule for vibrational excitation to be A v = 2, To 
the best of our knoujledge, ours is the first result on the CID 
process for the collinear (H, H 2 ) system making use of the 
accurate initio PES. A more detailed study is in progress. 
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APPENDIX A 


The Fourier transform F(tju) of a function f(x) is 

given by 

CXJ , 

F(uj) = / f(x) dx ..(A1) 

-GO 

uue evaluate the integral numerically, lue consider the 
discrete Fourier transform (OFT); 

F(u,) f(nAx) ..(A2) 

n=®0 

vuhere Ax is the sampling interval (mesh spacing). For equi- 
distant sampling points, Eq, (A2) becomes, 

F(»,A») - ix f(n^x) e-K^^^XnAx) _ 

n=®0 

117 

Standard softujares like the NAG depend only on the summa- 
tion indices and not the nature of the independent variaole 
Therefore lue rewrite Eq, (A3) in order for it to be compatible 

with the subroutine C06FCF as follows; 


^ ^ . -i2Ttinn/N 


F = Ax 2 f e 
n=0 ^ 

.. (AA) 

where Aw Ax = 2^/N; 

(A 5) 

F = F(mAw) and f_ = f(nAx) 
m n 

. . (A6) 



In tmo dimensions Eq, (A4) becomes 


F = Ax 2 g-ini>|n^2VN 

n-)=D 


N-1 


Ay 2 f ( n-| , 02 ) e 


n2»0 


-im2n227T/N 

*.(A7) 


We have to evaluate the DPT s along the columns and with this 
result, subsequently evaluate the OFT s along the rows. This 
approach is called the row.«column method because it can be 
viewed as equivalent to organising the input data into sets 
of row and column vectors in an array and performing the 
computation. 

If the values of the function are known from the point 
-N/2 to the point n/ 2 - 1 at equidistant points, Eq, (A3) becomes 

= Ax S ■ f(nAx) ..(A8) 

n- - n/ 2 
where -N/2 4m4^/2-1 
Change of indices 

k = m + n/2 5 1 = n +N/2 ..(A9) 

leads to 

N-1 ■ , -iCk-N/2)(l-N/2)2*VN 

F[(k-N/2)Aw] = Ax 2 f[(l-N/2)Ax] e 

1=0 


..(A10) 



which on simplification yieldsj 


F[(k-N/2)^uj] = gi(k-N/2)Tt f[(i.N/2)Ax]« 

1=0 

-i27r)<l/N 


..(A11) 


The second derivative of the function f(x) can be 
obtained using the properties of Fourier transform. The inverse 
Fourier transform of Eg. (A1) is 


00 


f(x) = -^ / F(uj) duj 

-00 

The second derivative of f(x) with respect to x is 


..(A12) 


*1 00 ■’ 
f”(x) = / -uj^F(u)) 8^ ^ duj 


27T 


..(A13) 


-00 


Discretising Fq. (13) we obtain 


/ ,^1 / *2 / V i(mi&w)x 

f*‘(x) ~ S -(mAuj) F(mAw)e 

m=0 


AUJ ..(A143 


J F(mAuj) ..(A15) 

• m=0 


3 N-1 

f"(nAx) = S m^ F(mAw) e 


ziAAl 'r J rUAu,^ 


m=0 



li 


..(A17) 


f 


n 



N-1 

2 

m=0 


2 

m 


F 

m 


„imn2^/N 

8 


In order to compute the second derivative, first the 

set has to be calculated using Eq. (4) ujith a knowledge 

of the set ff„}. We have made use of the subroutine C06FCF of 
n 
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MAG to compute The result has to be multiplied by 

-m , where m is defined as in Eq, (A5) . Then we have to perform 

an inverse Fourier transform* For this, ujb have made use of the 

1 17 

subroutine C06GCF of NAG in conjunction with the subroutine 
C06FCF. The result when multiplied by the constants as in 
Eq, (A17) , gives the desired second derivative. 

In tujo dimensions Eq, (A17) becomes 


f"(niAx, nAy^ 



N-1 

E 

k=0 


ji(kAWj^)(inAx) 


3 

-] 


N-1 2 2 

X 2 (k 1 ) F. 1 e 
1=0 


i( lAw ) ( nAy) 

y ..(Ala) 
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We have predicted the absorption spectrum of the transition state configurations for the collincar reaction H + H 2 
H 2 + H by solving the time-dependent Schrodinger equation for the nuclear motion. The spectral features have been inter- 
preted in terms of the build up of the probability density function in the inner repulsive wall of the potential-energy surface. 


1. Introduction 

It has been recognized in recent years [ 1 ] that it is 
possible to obtain the spectrum of the transition state 
(TS) ^ for elementary chemical reactions. Arrowsmith 
et al. [2,3] reported the emission spectrum of the TS 
for the reaction 

Na 2 + F->NaF + Na* (Rl) 

by monitoring the “wings'” to the sodium D-line emis- 
sion. Polanyi and Wolf [4,5] interpreted the wings as 
arising from the emission from Na* while the separat- 
ing NaF is still in the vicinity. Brooks and co-workers 
[ 6 - 8 ] have reported evidence for absorption from the 
intermediate configurations in the reactions 

K + HgBr 2 KBr + HgBr , (R2) 

K + NaCl KCl + Na . (R3) 

Such TS spectra can provide information on the colli- 
sion dynamics complementary to what can be learned 
from state-to-state studies [9]. 

It is of fundamental interest to predict and observe 

^ Permanent address: School of Studies in Physics, Vikram 
University, Ujjain, Madhya Pradesh, India. 

2 In partial fulfillment of the requirements for the degree of 
Doctor of Philosophy. 

We use the term transition state in a general sense to refer 
to all configurations intermediate between reactants and 
products. 


the wings (if any) to the Lyman-a absorption (emis- 
sion) line in H + H 2 collisions. Polanyi and co-workers 
[ 10 , 1 1 ] have computed the absorption spectrum for 
the TS for the coUinear reaction 

H + H 2 ^ H 2 + H (R4) 

by determining the probability density distribution 
for the TS on the chemically accurate Siegbahn-Liu- 
Tmhlar-Horowitz (SLTH) [12,13] ab initio poten- 
tial-energy surface (PES) using the quasiclassical tra- 
jectory (QCT) technique [14]. They assumed vertical 
transitions from the ground state to one of the low- 
lying excited states ^0 in geometry that corre- 
lates with H*(2p) + H 2 . A correlation diagram of dif- 
ferent states for H 3 is presented in fig, 1 of ref. [15]. 
They used model PESs based on the available ab initio 
data for the upper PES. They studied the effect of in- 
creasing relative translational energy of reactants and 
reagent vibrational excitation on the features of the 
absorption spectrum. As the collision energy increased, 
the intensity of the spectrum diminished and there was 
a considerable red-shift for the wings to the Lyman-a 
line. Subsequently, they extended [15] their study to 
collisions in three dimensions (3D) using a diatomics- 
in-molecules [16] (DIM) surface for the upper state 
and concluded that the qualitative features of the spec- 
trum remained while going from coUinear geometry 
to 3D. 

Recently, Engel et al. [17] have carried out time- 
independent quantum mechanical (QM) calculations 
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and predicted the absorption spectrum for the collin- 
ear reaction (R4) for two different vibrationai states 
(v - 0 and I ) of the hydrogen molecule and at colli- 
sion energies below the reaction threshold (0.22 eV) 
using the same ground- and upper-state surfaces as 
Polanyi and co-workers did [ 1 1 ]. The QM and the QCT 
results agree with each other in the basic features and 
trends in the spectrum. 

In this Letter, we report the results of a time- 
dependent wave meclianical calculation of the absorp- 
tion spectrurn for the same collinear reaction (R4) for 
u = 0 and I. at two different collision energies, 0.41 
and I.OOeV, that are well above the reaction threshold, 
thus complementary to the results of ref. [17]. 

2. Computational method 

The two-dimensional time-dependent Schrodinger 
equation 

[(-nyiUiXd^ldqj + 02/3^2)+ F(^,, ^2)]'!^ 

= ifi d^/dt , (I) 

where and q 2 are the scaled and skewed coordi- 
nates [18], has been solved using the explicit integra- 
tion method described elsewhere [19]. An L-shaped 
box in the space, 1.7 < 8.19 au and 0.4 

<^^2 ^ enclosing the interaction region has 

been divided into a larger number of mesh points. The 
values of the time step At (= 0.53875 X s) and 
the mesh spacing Ay ( j = ^ = 0. 1 1 au) have been 

chosen such that they meet the stability criterion for 
the explicit integration method [19]. The initial wave- 
function 'T(r = 0) was chosen as 

= 0) = Mcii. q.J^ <K^1 ) X^iq ,) , (2) 

where is the /hh state Morse oscillator wavefunc- 
tion for the hydrogen molecule. (I>(q ^ ) is the minimum 
uncertainty Gaussian wavepacket representing the rela- 
tive translational motion, given by 

<^(^/,) = (27r52)-l/4 

Xexp[-(^] - ■ O) 

The momentum space distribution is given by its 
Fourier transform: 


/’(*)= (26 2/,r)'/'> 

X exp[— (A: — A'q)*62 + iq^Q(k — Atq)] . (4) 

The average kinetic energy of this packet is 

(E> = (fi^/2(i)(kl + m^-). (5) 

We have chosen 5 = 0.25 au;( 7 jQ = 5,50 au; and Icq = 
5.75 and 9.27 au for (E) = 0.41 and 1.00 eV, respec- 
tively. The time evolution of the system on the SLTH 
surface described by eq.(] ) was terminated at the time 
when the plot of the reaction probability as a func- 
tion of time levelled off, as shown in fig. I. We have 
checked the accuracy of our calculation by checking 
the normalization condition for ^ throughout the 
time evolution of the system. In addition, we find that 
aif average of the time-independent collinear reaction 
probability [20] over the translational energy distribu- 
tion used for the wavepacket is in agreement with our 
reaction probability value for u = 0 at {£) = 0.41 and 
1.00 eV to within 10%. 

The value of d^j dq 2 over each area element 
represents the probability of the TS in that area ele- 
ment. The cumulative value of dq^ d^2 
the entire time evolution period for each area element 
represents the intensity of the transition at an energy 
equal to the energy difference between the ground- 
and the excited-state surfaces in that area element. 



Fig. l. Reaction probability as a function of time with each 
time step = 0.53875 X 10“^^ s, for a specific initial condition. 
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We have used the DIM surface described in ref. [15] 
for the upper surface. We have assumed a constant 
electronic transition moment as was done by the ear- 
lier workers f 1 1, 15,1 7 1. We have used a histogram 
width of 0.2 eV in reporting the spectrum. 

3, Results and discussion 

The resulting absorption spectra for H3 for (E) = 
0.41 and 1.00 eV foru = Oand I are reported in fig. 2. It 
is clear that there are substantial wings to the Lyman-o 
line at the lower energy region. Two prominent peaks 
have been identified as a and b. With an increase in 
collision energy, the intensity of the peak close to the 


line decreases and that of the peak in the lower- 
frequency region increases. Similar trends have also 
been reported by the earlier workers [1 1,15,17]. 

In fig. 3, the cumulative value of the probability 
density function (PDF) over the entire grid is super- 
imposed over the difference potential contours (AV 
= Fjj - the difference in the potential energy be- 
tween the ground- and the excited-state surfaces). 
Hence each configuration can be read as a spectral 
frequency corresponding to the energy difference 
AV. The greater the PDF at any particular AF, the 
greater the absorption intensity at the corresponding 
V. Thus the peaks in the absorption spectra can be 
identified with the accumulation of the PDF at a par- 
ticular range of AF values. The peak labelled “a” (at 


AV (eV) 



X>/I000.cm"’ 


Fig, 2. Absorption spectra for Hf under different initial consitions. Average collision energies are indicated against each spectrum. 
L^ represents the Lyman-o: line. 
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Fig. 3. Cumulative probability density function in the configuration space is indicated by dashed lines, along with numbers indicating 
its relative magnitude in different regions. The solid lines represent AC contours with the numbers indicating the values in eV. a 
and b indicate the PDF corresponding to the peaks in fig. 2. X indicates the saddle point. 


72600 env"” ^ ) in fig. 2 for (E) = 0.41 eV arises due to 
the probability density contour (PDC) marked 3 and 

5 running parallel to the r/ j axis in the AV region be- 
tween 8 and 9 eV. The peak labelled “b'' in the lower 
wavenumber region of the spectrum (at 5 1600 cm 
arises due to the build up of the PDF marked 3 near 
the saddle point. As the collision energy is increased 
to 1 eV, the PDC marked 5 is found to be missing in 
the AF region between 8 and 9 eV and hence the in- 
tensity of the peak “a" falls off. But there is a build up 
of tlie PDF marked 4 in tlie AF region between 5 and 

6 eV explaining the increase in the intensity of the 
peak ''b'’ with increase in (E). 

It is interesting to note that increase in v from 0 to 1 
does not affect the shape of the wings to any significant 
extent as was found by the earlier workers [1 1,15,17], 


4. Conclusion 

We have predicted the absorption spectrum for the 
TS for the collinear reaction {R4) using the SLTH sur- 


face for the ground state and a DIM surface for one 
of the upper states that correlates with H*(2p) + H 2 , 
using a constant electronic transition moment. Two 
prominent peaks have been identified and interpreted 
in terms of the build up of the PDF in the configura- 
tion space. The effect of including non-collinear colli- 
sions, a variable transition moment, and transition to 
the other low-lying upper electronic states and its im- 
plications to an experimental observation of the wings 
to the Lyman-o line are currently under investigation. 
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TIME-DEPENDENT QUANTUM MECHANICAL APPROACH TO 
REACTIVE SCATTERING 

V. MOHAN* and N. SATHYAMURTHY 
Department of Chemistry, Indian Institute of Technology, Kanpur 208016, India. 

ABSTRACT 

There has been a renewed interest in the time-dependent quantum mechanical 
approach to reactive scattering in the last few years. The different computational methods 
and applications that have created such an interest are reviewed here. 


1 INTRODUCTION 

* LTHOUGH it was realised more than fifty years 
ago that, in principle, any chemical problem 
could be solved by an appropriate formulation of 
it in quantum-mechanical terms*, in practice, 
developments in that direction had to await the 
arrival of modern electronic computers. In the 
early seventies one could calculate reaction prob- 
abilities ( F* 1 readily for exchange reactions of the 
type 

A + BC-*AB + C (Rl) 

in collinear geometries by solving the time- 
independent Schrodinger equation for the nuc- 
lear motion, given the electronic energy (for 
example, sec ref. 2). That exact (converged) sol- 
utions to the problem in three dimensions (3D) 
also could be obtained was demonstrated in 1975 
when Kuppermann and Schtaz’^^and Elkowitz 
and Wyatt’' published their results for the 
reaction 

H + (R2) 

at low energies, making full use of the symmetry 
of the system. There has not been any report on 
results of similar accuracy for any other reactive 
system in the past ten years. Even for the reaction 
(R2), some converged results have become 
available^ only recently on the chemically ac- 
curate Siegbahn-Liu-Truhlar-Horowitz’ (slth) 
potential-energy surface (pes) for the ground and 


• In partial fulfilment of the requirements for the degree of 

Doctor of Philosophy. 


the first excited vibrational states (v = 0.1) of the 
reactant. While some progress can be expected in 
the next few years, thanks to the supercom- 
puters/parallel processors becoming available, 
yet the enormity of the problem remains. In a 
parallel development, Mazur and Rubin** 
showed that the time-dependent Schrodinger 
equation (tdsi^ could be solved for a reaction of 
the type (Rl) using a model potential based on 
Harmuth integration scheme*. McCullough and 
Wyatt"* (mw) did an elaborate study in which they 
solved the tdse for the reaction (R2) in collinear 
geometries using a realistic pes®. Unfortunately ; 
their approach could not be used readily for 
many systems as (i) the computer memory and 
time requirements were large (ii) the resolution of 
the state-to-state F* information was poor and 
(iii) there was very little hope of the method being 
extended to three dimensions. In the last few 
years, however, the situation has changed drama- 
tically. Newer methods’-'® of solving the tdse 
have become available and they are an order of 
magnitude faster than that one used by mw. Thus 
there is a renewed hope of studying 3D atom- 
diatom exchange reactions quantum mechani- 
cally. Also Kulander" and Leforestier’^ have 
shown that it is possible to deconvolute the 
average reaction probabilities ( <F*> ) that result 
from a time -dependent quantum mechanical 
(tdqm) study. They have also shown that the 
tdqm approach can be used successfully in 
studying collision-induced dissociation (cid) 
processes -an area where the time-independent 
approach has been carried out successfully only 
for a model system’^ because of the continuumj 
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problem. There are also other attractions offered 
by the time-dependent approach. It is analogous 
to classical trajectories in that a pictorial rep- 
resentation of the chemical reaction is possible — 
giving us a “feeF’ for the molecular collision. It 
has been shown’ that flux patterns could be 
obtained for a chemical reaction on a pes and the 
whirlpools in the resulting picture can provide 
insight into reactive scattering resonances. Very 
recently we have also shown ^ ^ that by solving the 
TDSE, we can construct the probability density 
function and predict the absorption spectrum for 
the transition state of the collinear reaction (R2). 
Gas-surface scattering problems have also been 
shown to be amenable to solution via the tdse. 
Because of the renewed interest in the area, we 
present here a review of the existing methodology 
and applications of the tdqm approach to a study 
of reactive scattering processes. 


2. METHODOLOGY 
Mazur and Rubin^ solved the tdse 


ifT = rh 


5T 

dt 


by separating the real and imaginary parts in'F. 
Letting T = u -f iv, the equation becomes 



Hv = 



They approximated the spatial derivative in the 
hamiltonian and the temporal derivative by 
second order and first order finite differences 
respectively. Subsequent workers could deal with 
the complex waveftmction directly. Kulander^^ 
used a five-point finite difference formula for the 
spatial derivative and a fifth order predictor- 
corrector method for propagating the wavefunc- 
tion in time. Kosloff and Kosloff^^ have shown 
recently that the second order spatial derivative 
can be evaluated by a fast fourier transform (fft) 
of the wavefunction accurately. They used a 


second order finite difference formula for 
evaluating d^/dt. 

The time-evolution of the wavefunction can be 
expressed elegantly in terms of a time evolution 
operator U as 

T(f) = {/T(to) 

with 1/ = exp[ -iH{r-ro)/fi]- For practical 
reasons, this exponential operator has to be 
truncated. It was represented by mw as 

(/ = [1 ~ iH (f ~ to)/fi]/[l + iH (r ro)/h] 

which is unitary as is the untruncated U. But 
using it for determining the time evolution of a 
chemical system, involves inverting a very large 
matrix. Therefore the method is computationally 
expensive and it requires a large computer 
memory. As a result it has been applied to a few 
systems only^^ *^. 

Askar and Cakmak^ have proposed an explicit 
method in which U was truncated after the linear 
term. By combining the forward and the back- 
ward evolution operators, they could devise an 
algorithm which was an order of magnitude 
faster than the implicit method of mw. 

The first order explicit-method proposed by 
Askar and Cakmak^ is equivalent to the second 
order finite difference in time used by Kosloff and 
Kosloff^ It also turns out that the first order 
explicit is identical to a second order explicit 
method. (As a matter of fact, truncation upto an 
odd order in time is equivalent to truncation upto 
the next higher order.) Therefore, the errors are 
of third order. The stability criteria to be met by 
the choice of steps in time and space for the 
method have been worked out elsewhere 

The above methods have been applied to 
different collinear triatomic systems. The only 
exception is the study of Jackson and Wyatt^^ 
who studied the time evolution of a planar H 3 
system. The speed of the explicit method gives us 
the hope of extending it to atom-diatom ex- 
change reactions in three dimensions in the very 
near future. 

In solving the tdse for an atom-diatom ex- 
change reaction, the initial wavefunction is rep- 
resented by a product of the eigen-function for 
the vibrational state of the diatom (say, a 
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Morse oscillator) and a translational wavepacket 
which is gaussian in the momentum (fc) space, 
centered at ko corresponding to the classical 
relative translational energy of the reactants. The 
width (5) of the gaussian is chosen to be ‘optimaP 
( = 0.25 a.u. in our studies). Such a wavefunction, 
when time evolved, would yield an average 
probability <F^> for reaction into all the 
available product vibrational (r') states. If we are 
interested in a state (t;)- selected rate coefficient 
at a translational temperature T, the wavepacket 
can be chosen accordingly If an overall rate 
constant is desired, an appropriate summation 
over the v states can also be incorporated in the 
initial wavefunction. Recently, it has been shown 
that the > can be deconvoluted to obtain the 
dependence of on the translational energy 
^tam product state specified reaction 

probability 

3. APPLICATIONS 

3.1 Atom-diatom exchange reactions 

By far the predominant application of the 
TDQM approach has been to calculate <P^> for a 
variety of atom-diatom exchange reactions. 
MW*^ computed the <P^> at different energies 
for the reaction (R2) on the Porter-Karplus® pes. 
Agrawal and Raff used the explicit method to 
compute rate coefficients for the reaction 

T-^HD-^TH^D {R3) 

on a London-Eyring-Polanyi-Sato (leps) surface. 
The results were shown to be in good agreement 
with the time-independent results of Garrett et 
Values of <P^> computed^^ on the slth 
surface also have been shown to be in accord with 
the time-independent results^®. 

Kellerhals et aP® computed <P^> for the 
reaction (Rl) on a model leps surface and studied 
its dependence on the location of the saddle 
point. It became clear that <P^) was sensitive to 
changes in the pes. Stroud et aP^ investigated the 
influence of v on <P^ > and its sensitivity to 
features of the pes for the reaction 

He + Hi^HeH^+H (R4) 


using a diatomics-in-molecules and a slightly 
different spline-fitted ab initio surface. Studies by 
Zuhrt et al^^ using three different pess for the 
reaction 

(R5) 

also showed that <P^> was strongly dependent 
on the topological features of the pes. 
Sathyamurthy et aP^ were able to demonstrate 
the utility of cubic splines to interpolation of ab 
initio pess by showing that results on a splinefit- 
ted surface were identical with those obtained on 
the original analytic surface. Thareja and 
Satyamurthy^^ have recently used the tdqm 
approach in testing the quality of an analytic fit to 
the ab initio pes for the reaction (R2). When 
different fits were obtained for the same original 
ab initio surface, the fit that led to the same <P* > 
as on the slth was taken to be the best. 

In addition to yielding P* at different levels of 
resolution, the time-evolution method provides a 
valuable insight into the detailed nature of the 
collision dynamics. Typical 3D perspective plots 
of I'Pj^ in configuration space at different time 
intervals give a pictorial representation of the 
reactants transforming into products and the 
resulting picture can be compared directly with a 
plot of a family of trajectories “Structures” in 
these plots reveal the complexity of the dynamics 
and are indicative of indirect collisions. Using the 
FFT algorithm, Kosloff and KoslofP^ investi- 
gated the conditions for complex formation in 
-\-H 2 and its isotopic analogs for which the 
pes has a deep potential well in it. 

From a knowledge of ^(r, r), quantal flux 
patterns could be drawn, in analogy with the fluid 
flow in classical hydrodynamics. Plots of the 
quantal probability flux vector 

j(r,t)--Im['P*(r.t)V^(r,t)] 

P 

were made by MW"' and quantum whirlpools 
were identified. These have been investigated 
subsequently by Hirschfelder** who described 
them as quantized vortices around wavefunction 
nodes. 
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3.2 CoUishn-induced dmociation processes 

In contrast to the formidable problems faced 
by the time-independent approach in studying 
ao processes, the time-dependent approach is 
applicable to the cm problem with as much ease 
as for the exchange reaction. Ford et re- 
ported the first TDQM results for the reaction of 
the type 

A + + C (R6) 

using a model potential. Subsequently 
Kulander^ ^ reported the results for the reaction 

H + + H (R7) 

on a LEPS surface. Leforestier et 0/^2.27,28 
studied a variety of model cm reactions over a 
wide range of energies. By deconvoluting 
>to obtain the energy dependence of for 
a cm process, they have been able to obtain 
information on the threshold behaviour, vibra- 
tional enhancement and inhibition and their 
dependence on They have shed further light 
on the dynamical nature of cm processes by plots 
of I'PP and j at different time intervals. 

3.3 Spectroscopy of the transition state 

Recently, we^^ have made use of the tdqm 
approach to predict the wings to the Lyman-a 
absorption line in H+H 2 collisions by obtaining 
the time averaged j'Fp values and relating them 
to absorption intensity. Work is in progress in 
using this approach in predicting the emission 
spectrum for the system and spectra of 
transition states for a few other systems. The 
approach thus promises to be an effective tool in 
the emerging area of ‘spectroscopy of the trans- 
ition state’. 

3.4. Gas-surface scattering 

Atom-surface scattering plays a vital role in the 
investigation of solid surface properties. The 
rapid advancement in the recent past on the 
experimental side of the atom-surface encounters 
pose a challenge to theoretical models. Time- 
independent quantal calculations have not 
becomie practicable yet. Classical and semi- 


classical methods have been used extensively in 
simulating the experimental results. But they are 
of limited validity. For the first time Agrawal and 
Raff^^ showed that the time-dependent approach 
could be used to study gas-surface scattering by 
representing an incident atom by a wavepacket 
and using the explicit integration method for the 
time evolution of the system. Yinnon and 
KoslofP^ have been able to compute scattering 
intensities, resonance strengths and other attri- 
butes using the FFT algorithm. Gerber et aP^ have 
used the same approach to study the scattering of 
He atom from a Cu surface with isolated Ar 
impurities. Kosloff and Cerjan^^ have studied the 
desorption and scattering of a He atom from W 
and Pt surfaces using the same method. 

4. CONCLUDING REMARKS 

The TDQM approach to reactive scattering has 
been receiving considerable attention, of late. It 
scores over the time-independent approach in 
that there is no need for construction of channel 
hamiltonians and matching of the wavefunctions. 
There is no increase in the complexity of the 
problem or the computational time requirement 
with an increase in the total energy of the system. 
As a matter of fact, there is a reduction in the 
computer time with an increase in energy- 
analogous to the classical approach. In contrast 
to the time-independent approach in which 
has to be computed for different and v and 
averaged over the appropriate distributions in 
order to calculate state-selected and overall rate 
constants, the tdqm method could be made to 
yield them directly by a suitable choice of initial 
T. Also, the deconvolution methods have made it 
possible to obtain in a single calculation P^ over a 
wide range of E^^- 

The time is not far off when the time- 
dependent approach could be used with much 
more ease than its time-independent analog in 
predicting rate observables for reactions in three 
dimensions. 
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